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CONVERSATIONAL HIGH-RESOLUTION MASS SPECTROGRAPHIC 


DATA REDUCTION 

Michael P. Romiez 
Ames Research Center 


SUMMARY 


A FORTRAN IV program is described which reduces the data obtained from a high-resolution 
mass spectrograph. The program (1) calculates an accurate mass for each line on the photoplate, and 
(2) assigns elemental compositions to each accurate mass. The program is intended for use in a 
time-shared computing environment and makes use of the conversational aspects of time-sharing 
operating systems. 


INTRODUCTION 


Since the organic chemist first realized the potential of high-resolution mass spectrometry as 
an analytical tool (ref. 1), a great effort has been spent in developing computer software to process 
the large amounts of data generated (refs. 2, 3, and 4). The first programs were written for use on 
large batch computer systems. Recently, however, the emphasis has shifted to programing for 
relatively small, dedicated computer systems, with the data acquisition and reduction in real-time 
(refs. 5 and 6). 

Obviously, a real-time data system has many advantages, the best of which is immediate data 
analysis and presentation. In the off-line schemes, typically a day or more separated the data 
acquisition and the computer analysis of the mass spectra. A disadvantage of the real-time systems is 
that it takes approximately 20 seconds (ref. 6) (for resolution of one part in 10,000) to acquire a 
mass spectrum. If greater resolution than one part in 10,000 is desired, even more time is neccessary 
to obtain the spectrum. For certain applications, for example, a directly coupled gas-chromatograph 
mass spectrometer, the sample concentration may be so small that it does not remain in the mass 
spectrometer long enough for the complete mass spectrum to be acquired. For these kinds of 
applications a photodetection system, which is an integrating detector, is superior to the real-time 
system. 

In the mass spectrometers used in real-time systems (ref. 7), the beam of ions comes to focus 
at a point. The magnetic field or ion accelerating voltage of the mass spectrometer must be scanned 
to sweep all the ion types past an electron multiplier at the point of focus. For samples of small 
concentration, as mentioned above, there may not be enough time to complete the scanning 
process. On the other hand, in a mass spectrograph (ref. 7), the ion optics of the instrument bring 
the sample beam to focus in a plane, and a photoplate placed in this plane can continuously record 
all the ions of the entire mass spectrum. With this scheme, even for small sample quantities, a 
complete mass spectrum can be obtained. For most of the analytical problems encountered in this 
laboratory, a photoplate recording system is the preferred method of data acquisition. 


The computer program described below is designed to reduce data acquired off-line from a 
photoplate, but, because it runs on a time-sharing computer, it has some of the advantages of a 
real-time program. A day or more still separates acquisition of the mass spectrum and its computer 
analysis, but after the program is in use, it communicates with the user at a conversational terminal 
and accepts various inputs from the user during execution. As will be discussed later, parameters 
that affect the calculations in the program must be selected on the basis of various characteristics of 
each photoplate. With this conversational program, the selections can be made and changed while 
the program is running, the same as in a real-time system. 


SYSTEM OVERVIEW 


The high-resolution mass spectrograph (CEC 21—1 10B) employed at Ames Research Center is 
of. the Mattuch-Herzog geometry. The sample to be analyzed is introduced into the ion source 
chamber of the mass spectrograph, and is ionized and fragmented by a collimated beam of ther- 
mally emitted 70 volt electrons. The positively charged fragments of the sample are accelerated 
through a focusing slit and passed into a magnetic field where ion-type groupings separate according 
to their mass. Each charged fragment type forms a slit image on a photoplate placed in the focal 
plane of the instrument. The photoplate can record up to 45 exposures (mass spectra) of typically 
several hundred slit images (peaks) each. 

Each mass spectrum is transcribed from the photoplate to seven track BCD magnetic tape by a 
custom-designed comparator-microphotometer interfaced with a tape writing unit. The exposed 
photographic plate is scanned past the optical density measuring photomultiplier on a precision 
driven table. When the circuitry senses a line on the photoplate, the table is slowed to tape writing 
speed, and optical transmittance data are passed to the digitizing circuitry every one-half micron. 
After the last optical density reading in a slit image profile is recorded, the distance of this last 
reading, in half microns from the low mass end of the photoplate, is encoded. 

At the conclusion of one complete spectrum transcription, the data tape will contain a profile 
and a distance value for each peak in the mass spectrum. This mass spectrograph data tape is read 
into disk storage on the Ames Research Center, IBM 360 Model 67 computer under an IBM sup- 
plied time-sharing operating system (TSS) and becomes the input file for the data reduction 
program. 

From the shape of the line profiles and the distances of the last optical density readings in each 
profile, the position of the center of each peak can be determined to within an error of less than 
one-half micron. The masses of each ion type are related to each other by the square of the distance 
that separates them. Peaks of known mass are interspersed in the sample spectra by having a mass 
standard (perfluoroalkane) always present in the ionizing chamber. These peaks of known mass are 
used to calculate a distance to mass relationship for the length of the photoplate. From this 
relationship, the masses of the ion types present in the sample spectrum can be determined to 
within an error of a few parts per million. This mass measurement accuracy is, for the vast majority 
of problems, sufficient to assign unambiguous elemental compositions to each fragment in the 
sample spectrum. The final output is a tabular listing of the accurate mass and elemental composi- 
tion of each peak in the mass spectrum. 
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COMPUTER PROGRAM DESCRIPTION 


The computer program can be logically divided into three sections: peak center calculation, 
standard mass identification, and accurate mass calculation and element fitting. 


Peak Center Calculation 

Because the accurate mass calculated for each peak in the mass spectrum can only be deter- 
mined within the error of the peak center calculations, the determination of the center of each peak 
profile is the most important aspect of the data reduction process. For an ideal line profile, whose 
shape is a gaussian distribution curve, the accurate center is easily discerned. In general, however, 
not all profiles in a spectrum will have a gaussian shape. There will be profiles with complicated 
shapes that describe lines which are not completely separated, and there will be wide profiles with 
flat tops that correspond to overexposed lines on the photographic plate. In addition, always 
superimposed on each profile will be a certain amount of random error commonly referred to as 
noise. An algorithm for peak center calculation must first attempt to filter out mathematically as 
much random error as possible, without degrading the underlying profile data, before addressing 
itself to the determination of the center of each peak. 

Data smoothing— On examining a plot of a peak profile, one would tend to draw through the 
points, a line which best fits them, thereby obtaining the best smooth curve through the points. 
Numerically, this can be done by fitting a least squares polynomial to the data profile. 

The set of points to be smoothed is fit to some polynomial curve. 

N 

P(x) = ^ ajxi 
j=0 

The coefficients designated a are computed so that when the x values are substituted into the 
polynomial equation, the square of the differences between the computed values and the profile 
data points are a minimum. Methods used to calculate the a coeffi- 
cients are described in most textbooks on numerical analysis (refs. 8 — 

and 9). The degree of the polynomial (which can be varied from 1 to 
10) and the number of points used (2 to 200) in the moving polyno- 
mial smoothing function can be selected and changed while the pro- 
gram is executing. The sketch shows how the smoothing process 
works, for an overlapping 1 1 point polynomial. 

As in the first grouping in the sketch a new value for the central 
position in the group is calculated from the least squares polynomial. 

The point at the left is dropped and a new point added from the 
right. Again a central point is calculated from the new polynomial 
and the process is repeated until a smoothed curve is obtained for the 
entire profile. An intermediate grouping is also indicated in the 
sketch. All peak profiles in the mass spectrum are smoothed in this 
fashion. 
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Peak center calculation— The first derivative of each peak profile is calculated using over- 
lapping polynomials as above. From the shape of the first derivative curve, the number of peaks and 
the center of each peak in the profile can be determined (ref. 1 0). 


Not all line profdes are analyzed in this fashion. Narrow profiles which can only contain a 
single peak are analyzed by a scheme that uses less computer time. If the total number of points in 
the peak profile is below some minimum number (selected during the program execution), the 
profile is best fit to a gaussian curve and the center and height determined directly from this curve 
(ref. 4). 
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This second degree polynomial can be fit to the entire profile to yield directly the center and height 
of a peak. 


Occasionally, a peak is so overexposed on the photoplate that its shape cannot be properly 
analyzed by the methods described above. For this contingency a provision exists whereby a profile 
with an excessively large number of data points (the number is selected during program execution) 
will be displayed for editing on the conversational terminal. Long leading edges and tails (which are 
• a phenomenon caused by the emulsion on the photoplate during an overexposure) can be elimi- 
nated and the profile can otherwise be edited before any further analysis. 

Peak intensity calculation — With electrical ion detection, the abundance of each ion type can 
be accurately determined, but the intensity measurement capabilities of a photoplate are very 
limited (ref. 1 1). Because of this limitation, no attempt is made at accurate ion abundance measure- 
ments. The largest intensity reading in each profile is used to calculate an approximate relative 
abundance for each ion type. 

A list of peak centers, in millimeters from the left edge of the photoplate, and high intensity 
values for each peak in the mass spectrum are passed to the standard mass identification portion of 
the program. 
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Standard Mass Identification 


As mentioned above, not all the lines in a spectrum are due to the analyzed sample; some arise 
from an internal mass standard (perfluoroalkane). It is the task of this section of the program to 
identify which peaks correspond to masses from the mass standard and to use these to establish a 
distance-to-mass relationship for the length of the photoplate. 

Mass identification — To begin with, the distances of the first 40 peaks in the list are displayed 
on the conversational terminal. At this point the program must be supplied the known mass of at 
least two lines in the table. Using these two (or more) lines, the program will calculate distances 
where it expects to find 160 of the most prominent peaks in the perfluoroalkane spectrum. Each 
time the program is successful in its search, it incorporates the new line into its standard mass table. 
At the end of one search iteration, the program will use the table of identified masses to calculate a 
distance-to-mass relationship. 

Distance to mass conversion — Theoretically, there is a linear relationship between the square 
root of the mass of each ion type and the distance that separates the ions (ref. 7). In practice, 
because of inhomogeneous areas in the magnetic field used for separation, the relationship is not 
exactly linear, but is best approximated by overlapping polynomial functions that are fitted to the 
standard mass values spread throughout the length of the photoplate. The degree of and the number 
of points spanned by the polynomials may be selected during the program execution. Generally, the 
best results are obtained from a second degree polynomial fit to four standard masses. The coeffi- 
cients of the polynomial are calculated and saved for each overlapping set of calibration masses. To 
assess the accuracy of the conversion relationship, these polynomials are used to recalculate the 
masses of the identified standard mass lines, and this information is displayed at the conversational 
terminal. The difference between the calculated mass and the known mass of the lines of the 
calibration compound gives a very good idea of the accuracy of the distance-to-mass relation- 
ship in the area of interest. At this point misidentified peaks in the list may be deleted, and missed 
peaks may be added from the conversation terminal. If there has been any input from the user, the 
program loops back to the standard mass search routine and a new conversion relationship is 
calculated. This process may be repeated until one is sure that the calculated distance-to-mass 
relationship is correct. 

Element fitting parameters — If there are no additions to or deletions from the standard mass* 
list, the program prompts the user for element fitting parameters. A tolerance for element fitting 
must be chosen; it is the variance in millimass units between the calculated mass and the mass of an , 
acceptable elemental composition. Next the elements considered for possible compositions must be 
selected. Carbon, hydrogen, nitrogen, and oxygen are always included in the candidate element 
table, and as many as seven additional elements may be added for consideration at execution time. 

After the element fitting parameters are entered, the user may choose either (1) to write all 
information calculated thus far into a disk file, where it may be accessed by a program that uses the 
computer center line printer for output, or (2) to let the program continue with the final results 
displayed at the conversational terminal. Usually there will be too much information for listing on 
the relatively slow conversational terminal. 
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Mass Calculation and Element Fitting 


Mass calculation — The distance-to-mass relationship just obtained, is used to calculate an 
accurate mass for all the lines in the mass spectrum. Each mass is calculated using the polynomial 
that was fit to the nearest group of calibration masses. Next, all masses corresponding to ion species 
commonly found in the background of the mass spectrometer and all ion types of the marker 
compound are eliminated from the list of calculated masses. It is this abbreviated list of accurate 
masses that represents the mass spectrum of the unknown sample. 

Element fitting — All possible permutations of the elements submitted for element fitting are 
considered as compositions for each of the masses calculated. The minimum number of each isotope 
considered is zero. The maximum number varies as follows: for carbon-12 it is the integer resulting 
frofn the division of the accurate mass by 12; for carbon- 13 it is one; for hydrogen it is twice the 
maximum number of carbon atoms plus eight; for nitrogen-14 it is two; for oxygen-16 it is three. 
All, other isotopes are entered during program execution and the maximum number of each is one 
of the input parameters. All permutations whose masses fall within the millimass unit tolerance 
specified are displayed along with the calculated accurate masses. 

Output — The final output consists of a relisting of the standard mass table of part II and a 
table of all calculated masses and their associated elemental compositions. If no permutations of 
elements fall within the tolerance, the calculated mass is listed with an appropriate comment. There 
is no provision for special handling of multiply-charged ions in the mass spectrum. Generally, if an 
obviously multiply-charged species (calculated to be in the vicinity of half an integral mass unit) is 
doubled, the corresponding singly-charged ion can be found at this value. The accurate mass of the 
multiply charged ion species is listed with no associated elemental compositions. 


CONCLUSIONS 


A typical standard mass listing is shown in appendix A. The largest deviation from the exact 
mass in the entire list is less than five parts per million. This accuracy compares favorably with 
previously reported high resolution data analysis programs. In addition, for directly coupled gas- 
* chromatograph mass spectrometery applications, a photoplate detection system used with the con- 
versational data reduction system described here is more sensitive and more efficient than pre- 
' viously reported systems. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 94035, March 22, 1973. 
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APPENDIX A 


SAMPLE OF STANDARD MASS ACCURACY 



DISTANCE 

CALCULATED 

EXACT MASS 

DIFFERENCE 

DIFFERENCE 



MASS 


IN MMU 

IN PPM 

1 

15-3A97 

26.01569 

26.01565 

0.04 

1.A7 

2 

16.9022 

27.023A2 

27.023A7 

-0.06 

2.0A 

3 

18.3888 

28.00616 

28.00615 

0.01 

0.29 

A 

22.7609 

30.998AA 

30.998A0 

0.04 

I.3A 

5 

2A.1629 

31-98985 

31.98983 

0.02 

0.51 

6 

29.5927 

35.97659 

35.97668 

-0.09 

2.AA 

7 

3A.7268 

39.96239 

39-96239 

0.00 

0.02 

8 

39.6593 

A3. 98983 

A3. 98983 

-0.00 

0.02 

9 

A6.6IA7 

A9- 99686 

A9. 99680 

0.06 

1.10 

10 

A7.7399 

51 .OOA65 

51 . OOA63 

0.02 

0.A3 

11 

52.0936 

5A. 99833 

5A.998A0 

-0.07 

1.19 

12 

59.3566 

61.9968A 

61.99680 

0.04 

0.71 

13 

67.1731 

69.998A1 

69.99856 

-0.15 

2.20 

1A 

70.9115 

73.99688 

73.99680 

0.08 

1.10 

15 

77.218A 

80.99519 

80.99521 

-0.02 

0.25 

16 

81.5609 

85.99656 

85.99681 

-0.25 

2.90 

17 

87.A290 

92.995AA 

92.99521 

0.23 

2.52 

18 

93.0803 

99.99501 

99.99521 

-0.20 

2.00 

19 

96.9997 

10A. 99532 

10A.9952I 

0.11 

1 .OA 

20 

102.3302 

111.99335 

111.99361 

-0.26 

2.33 

21 

106.0392 

116.99518 

116.99521 

-0.03 

0.26 

22 

111.0970 

123.99365 

123.99361 

0.04 

0.3A 

23 

116.0157 

130.991A1 

130.99202 

-0.61 

A. 68 

2A 

116.7109 

131.99539 

131.99537 

0.02 

0. 18 

25 

122.82A0 

1A0. 99512 

1A0. 99521 

-0.09 

0.61 

26 

12A.15A6 

1A2.99197 

1A2.99201 

-0.04 

0.28 

27 

126.7720 

1A6. 96262 

1 A6.962A8 

0.14 

0.95 

28 

128.7AA3 

1 A9.99033 

1A9.990A1 

-0.08 

0.5A 

29 

131.9597 

15A. 99195 

15A. 99201 

-0.06 

0.A1 

30 

136.3722 

161 .99051 

161.990A1 

0.10 

0.62 

31 

139.A679 

166.99177 

166. 99201 

-0.24 

1.A6 

32 

1A0.6919 

168.98871 

168.98882 

-0.11 

0.62 

33 

1A3.72A6 

173.99051 

173.990A1 

0.10 

0.59 

3A 

1A6.7139 

178.9916A 

178.99201 

-0.37 

2.08 

35 

IA6.8967 

180.98895 

180.98882 

0.13 

0.70 

36 

150.8293 

185.99030 

185.990AI 

-0.11 

0.61 

37 

15A.8685 

192.98889 

192.98882 

0.07 

0.37 

38 

158.8352 

199.98709 

199.98722 

-0.13 

0.65 

39 

161 .6289 

20A.9889A 

20A. 98882 

0. 12 

0.58 

AO 

165.A817 

211.9872A 

21 1 .98722 

0.02 

0.08 

A) 

168.1976 

2 16. 98875 

216.98882 

-0.07 

0.30 

A2 

I71.9A56 

223.98732 

223.98722 

0. 10 

0.A5 

A3 

17A.5880 

228.98836 

228.98882 

-0.46 

2.01 

AA 

175.6358 

230.98565 

230.98562 

0.03 

0.1A 

A5 

178.2386 

235 - 9872A 

235.98722 

0.02 

0. 10 

A6 

180.815A 

2A0. 98862 

2A0. 98882 

-0.20 

0.82 

A7 

181.8369 

2A2. 98585 

2A2. 98562 

0.23 

0.95 

A8 

18A.376A 

2A7. 98717 

2A7. 98722 

-0.05 

0.22 

A9 

187.8883 

25A. 98563 

25A. 98562 

0.01 

0.0A 

50 

190.3679 

259.987A5 

259.98722 

0.23 

0.90 

51 

191.3517 

261.98396 

261 .98A02 

-0.06 

0.2A 

52 

193.8005 

266.98578 

266.98562 

0. 16 

0.59 

53 

I9A.7711 

268.98230 

268.982A2 

-0. 12 

0.A5 

5A 

197.1875 

273.98A13 

273.98A02 

0. 1 1 

0.A1 

55 

199.5817 

278.98529 

278.98562 

-0.33 

1.17 

56 

200.5325 

280.982A3 

280.982A2 

0.01 

0.0A 

57 

202.8978 

285.98399 

285.98A02 

-0.03 

0.11 

58 

205.2A27 

290.98575 

290.98562 

0.13 

0.A6 

59 

206.1733 

292.98258 

292.982A2 

0.16 

0.5A 

60 

208.A91A 

297.98A02 

297.98A02 

-0.00 

0.01 

61 

211.7035 

30A.982A1 

30A.982A2 

-0.01 

0.03 

62 

217.1272 

316.98269 

316.982A2 

0.27 

0.86 

63 

218.0198 

318.97910 

318.97923 

-0.13 

O.AI 

6A 

220.2AA2 

323.98096 

323.98083 

0.13 

0.39 

65 

222.A508 

328.9818A 

328.982A2 

-0.58 

1.76 

66 

223.3280 

330.97922 

330.97923 

-0.01 

0.03 

67 

225.51 1A 

335.98091 

335.98083 

0.08 

0.25 

68 

227.6800 

3A0.98186 

3A0.982A2 

-0.56 

1.65 

69 

228.5A19 

3A2. 97977 

3A2. 97923 

0.54 

1.59 

70 

230.6881 

3A7.980A7 

3A7. 98083 

-0.36 

1.03 

71 

233.6681 

35A. 97938 

35A. 97922 

0.16 

O.AA 

72 

235.7785 

359.98078 

359.98083 

-0.05 

0.13 

73 

238.7097 

366.97938 

366.97923 

0.15 

0.A2 

7A 

239.5A06 

368.97603 

368.97603 

0.00 

0.01 

75 

2A3.6715 

378.97910 

378.97923 

-0.13 

0.3A 

76 

2AA.A89A 

380.9759A 

380.97603 

-0.09 

0.2A 

77 

2A6.5286 

385.97791 

385.97763 

0.28 

0.72 

78 

2A9.3603 

392.97586 

392.97603 

-0.17 

0.A2 

79 

25A.160A 

AOA.976A5 

AOA. 97603 

0.42 

1.05 
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DISTANCE 

CALCULATED 

HASS 

EXACT MASS 

DIFFERENCE 
IN MMU 

DIFFERENCE 
IN PPM 

80 

257.3239 

412.98206 

412.98242 

-0.36 

0.86 

81 

258.8923 

416.97604 

416.97603 

0.01 

0.02 

82 

263.5553 

428.97626 

428 . 97603 

0.23 

0.54 

83 

264.3252 

430.97260 

430.97284 

-0.24 

0.55 

84 

266.2459 

435.97467 

435.97444 

0.23 

0.52 

85 

268.1541 

440.97582 

440.97603 

-0.21 

0.48 

86 

268.9139 

442.97233 

442.97284 

-0.51 

1.14 

87 

273.6419 

454.97284 

454.97284 

0.00 

0.01 

88 

277.9075 

466.97284 

466.97284 

0.00 

0.00 

89 

282.3171 

478.97289 

478.97284 

0.05 

0.10 

90 

283.0455 

480.96973 

480.96964 

0.09 

0.20 

91 

287.3917 

492.96960 

492.96964 

-0.04 

0.07 

92 

291.6849 

504.96941 

504.96964 

-0.23 

0.46 

93 

295.9264 

516.96986 

516.96964 

0.22 

0.43 

94 

300.1182 

528.96730 

528.96964 

-2.34 

4.43 

95 

300.8127 

530.96465 

530.96445 

0.20 

0.38 

96 

304.9507 

542.96604 

542.96645 

-0.41 

0.75 

97 

309.0428 

554.96634 

554.96645 

-0.11 

0.21 

98 

313.0897 

566.96659 

566.96645 

0.14 

0.25 

99 

317.0956 

578.96471 

578.96645 

-1.74 

3.01 

100 

317.7595 

580.96341 

580.96325 

0.16 

0.27 

101 

321.7162 

592.96277 

592.96325 

-0.48 

0.80 

102 

325.6343 

604.96374 

604.96325 

0.49 

0.80 

103 

329.5140 

616.96308 

616.96325 

-0.17 

0.27 
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APPENDIX B 


PROGRAM LISTING 


C JANUARY 9,1973 MICHAEL ROMIEZ 

C 

C THIS PROGRAM WILL CONTROL THE READING OF THE AMES HIGH RESOLUTION DATA TAPE INTO A 
C DISK FILE WHICH CAN BE ACCESSED BY THE MAIN DATA REDUCTION PROGRAM. 

C 

C THE DATA ON THE TAPE IS ENCODED IN SEVEN TRACK BCD FORMAT; IT IS CONVERTED INTO EBCDIC 
C BY AN AMES LIBRARY SUBROUTINE ALBYTES. THE DATA IS CONVERTED TO INTEGER FORMAT BY AN 

C ASSEMBLY SUBROUTINE TAPERD. THE RAW DATA IN INTGER FORMAT IS WRITTEN INTO A DISK FILE. 

C 

C VARIABLES 

C 

C I BYTES 

C 

C I HOLD (300) 

C 

C INTTAB(400) 

C PNTPER (40) 

C 

C D I STAB (AO) 

C PEKNUM 

C TOTP 

C ERR 

C 

C RAW I NT (25000) 

C NUMINT (1000) 

C 

C RAWDIS(IOOO) 

C NEXINT 

C NEXDIS 

C NOL I ST 

C 
C 

c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

LOGICAL NOLIST 

I NTEGER D I STAB , PNTPER , PEKNUM , TOTP , ERR 
INTEGER RAW I NT , NUM I NT , RAWD I S 
COMMON/AMESTP/ I BYTES , I HOLD (300) 

COMMON/MASPEC/ I NTTAB (AOO) , PNTPER (40) ,DISTAB(AO) .PEKNUM, TOTP, ERR 
COMMON/F I LE/RAW I NT (25000) , NUM I NT ( 1 000) , RAWD I S ( 1 000) 

C CLEAR ALL COMMON AREAS 
DO 500 N=1 ,25000 

500 RAW I NT (N)=0 

DO 501 N=1 ,1000 
NUM I NT (N)=0 

501 RAWD I S (N)=0 

C NOLIST IS THE LOGICAL VARIABLE TO OPT FOR A PRINTER LISTING OF THE 
C TAPE FILE 

READ (5,50) NOLIST 
NEX I NT=1 
NEXDIS-1 

1000 DO 1001 1-1,300 

1001 I HOLD ( I )-0 

C BYTES IS AN AMES LIBRARY SUBROUTINE THAT READS ONE RECORD OF A TAPE 
C FILE AT A TIME INTO THE I HOLD ARRAY IN ZONE DECIMAL FORMAT. I ERR IS 
C AN ERROR FLAG: 0- NORMAL RETURN; 1- END OF FILE; 2-1/0 ERROR 
CALL BYTES (I HOLD, I BYTES, I ERR) 

IF(IERR.EQ.I) GO TO 100 
I F( I ERR. EQ.2) GO TO 1000 

C TAPERD IS AN ASSEMBLY SUBROUTINE TO UNBLOCK THE DATA TAPE FORMAT 
C AND CONVERT THE ZONE DECIMAL NUMBERS TO INTEGER. ERR-0 FOR NORMAL 
C RETURN; ERR-1 FOR TAPE FORMAT ERROR 

1 CALL TAPERD 

I F (ERR. EQ. 1 ) GO TO 1000 

2 IF(. NOT. NOLIST) GO TO 300 

3 IF (PEKNUM. EQ.O) GO TO A 
LIMIT1-NEXINT+T0TP-1 

L I M I T2-NEXD I S+PEKNUM- 1 
K0UNT1-0 

DO A001 M-NEX INT,L IMIT1 
KOUNT1-KOUNT1+1 

4001 RAWINT (M)-INTTAB(KOUNTl ) 

NEXINT-LIMIT1+1 

K0UNT2-0 

DO A002 M-NEXDIS.LIMIT2 

K0UNT2-K0UNT2+1 

NUMINT (M)-PNTPER (K0UNT2) 

A 002 RAWD I S (M ) — D I STAB (K0UNT2) 

NEXD I S-LIMIT2+1 
GO TO 1000 

C UPON READING AN END OF FILE, THE CONTENTS OF THE RAW DATA COMMON 
C BLOCK ARE READ OUT INTO A DISK FILE NAMED IN THE DDEF STATEMENT 
A NEXINT-NEXINT-1 
NEXD I S-NEXD I S-l 


- THE NUMBER OF BCD CHARACTERS FOUND IN ONE TAPE RECORD BY SUBROUTINE 

- ALBYTES 

- THE STRING OF CHARACTERS FOUND BY ALBYTES PACKED FOUR TO A WORD IN 

- ZONE DECIMAL FORMAT 

- TABLE OF INTEGER INTENSITY VALUES RETURNED BY SUBROUTINE TAPERD 

- TABLE OF THE NUMBER OF INTENSITY READINGS PER PEAK IN THE INTTAB 

- LIST 

- TABLE OF DISTANCE VALUES. ONE VALUE FOR EACH PEAK IN PNTPER 

- TOTAL NUMBER OF PEAKS IN THE TAPE RECORD 

- TOTAL NUMBER OF INTENSITY READINGS IN THE TAPE RECORD 

- AN ERROR FLAG RETURNED BY SUBROUTINE TAPERD. ERR-0 NORMAL RETURN 

- ERR-1 FOR TAPE FORMAT ERROR 

- TABLE OF INTENSITY VALUES FOR THE ENTIRE TAPE FILE 

- TABLE OF THE NUMBER OF INTENSITY READINGS PER PEAK IN THE ENTIRE 

- TAPE FILE 

- TABLE OF DISTANCE READINGS. ONE FOR EACH PEAK IN NUMINT 

- POINTER TO NEXT VARIABLE IN THE RAWINT ARRAY 

- POINTER TO NEXT VARIABLE IN NUMINT AND RAWD IS ARRAYS 

- LOGICAL VARIABLE TO OPT FOR RAW DATA LISTING. TRUE FOR NO LISTING 

- AND FALSE FOR A LISTING 


C JANUARY 9,1973 MICHAEL ROMIEZ 

C 

WRITE (7) NEXINT,(RAWINT(l),l=l,NEXINT),NEXDIS,(NUMINT(l),RAWDIS(l),l=l,NEXDIS) 
STOP 

C IF AN END OF FILE HAS BEEN READ THE VARIABLE PEKNUM IS SET TO ZERO 
100 PEKNUM=0 
GO TO 2 

C THIS SECTION WILL PRINT THE CONTENTS OF THE AMES DATA TAPE IF 
C REQUESTED BY SETTING NOLIST EQUAL TO FALSE 
300 IF(PEKNUM.EQ.O) GO TO 4 
LENGTH* I BYTES/4+1 
1NTLN* I BYTES/3 
WRITE (6, 16) 

WRITE (6,10) I BYTES 

WR ITE (6, 1 1 ) (IHOLO (J) ,J=1 .LENGTH) 

WRITE(6, 17) 

WRITE (6, 12) (INTTAB(J) ,J=I , INTLN) 

WRITE (6,13) PNTPER 
WRITE (6, lit) DI STAB 
WRITE(6, 15) PEKNUM 
GO TO 3 

10 FORMAT (IX, 1 THE NUMBER OF CHARACTERS IN THIS RECORD ARE ',14) 

U FORMAT (1X.30A4) 

12 FORMAT (' INTENSITIES' ,3014) 

13 FORMAT (IX, 'PNTPER ',20 1 4) 

14 FORMAT (1 X , ' D I STAB ',1018) 

15 FORMAT(IX, 'PEKNUM', 14) 

16 FORMAT (I HO) 

17 FORMAT (1HX) 

50 FORMAT (L5) 

END 
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* JANUARY 9,1973 


MICHAEL ROM l EZ 


* THIS SUBROUTINE WILL CONVERT AN AMES HIGH RESOLUTION TAPE RECORD 

* IN EBCDIC FORMAT TO INTEGER FORMAT FOR THE MAIN DATA REDUCTION 

* PROGRAM TO USE. 

* 

* xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

TAPE CSECT 

ENTRY TAPERD 

USING TAPERD, 15 
TAPERD SAVE (14,12) 

LR 2,15 

DROP 15 

USING TAPERD, 2 

* LOAD AMES COMMON BLOCK ADDRESS INTO REG 3 

L 3.COM1 

* LOAD MASPEC COMMON BLOCK ADDRESS INTO REG 4 

L A, COM2 

* THIS LOOP WILL CLEAR MASPEC COMMON BLOCK 

LR 5 , 1 * 

L 6,-F'482' 

L 7,-F'O' 

CLEAR ST 7,0(5) 

A 5,«f'4' 

BCT 6, CLEAR 

* LOAD REGISTER 5 WITH I BYTES 

L 5, 0(0, 3) 

* ADVANCE AMESTP ADDRESS A FULLWORD TO BEGINNING OF I HOLD ARRAY 

A 3,«F'V 

* INITALIZE THE TOTAL BYTE COUNTER IN REGISTER 5- THIS IS USED TO TEST 

* FOR A NORMAL READ WITH THE FOLLOWING LOGIC. IF WE SUBTRACT SIX FROM 

* THE TOTAL BYTES INITALLY AND THREE EVERY TIME WE READ AN INTENSITY 

* OR SIX EACH TIME WE READ A DISTANCE THEN WHEN THE COUNTER EQUALS ZERO 

* THERE SHOULD ONLY BE SIX ALPHABETIC CHARACTERS LEFT IN THE ARRAY. 

* IF THIS IS NOT THE CASE WE TELL THE CALLING PROGRAM THAT THERE IS A 

* FORMAT ERROR 

S 5,-F'6' 

* LOAD REGISTER SIX AS THE TOTAL OPTICAL DENSITY READING COUNTER 

L 6,-F'O' 

* LOAD REGISTER 8 AS THE OP I CAL DENSITY PER PEAK COUNTER 

L 8,-F'O' 

* LOAD REGISTER 9 WITH THE ADDRESS OF PNTPER 

LR 9,4 
A 9, “F' 1600' 

* LOAD REG 10 WITH THE ADDRESS OF D I STAB 

LR 10,9 

A 10,-F' 160' 

* LOAD REG 11 WITH THE ADDRESS OF PEKNUM 

LR 11,10 

A 1 1 ,-F ' 1 60* 

* INITALIZE REG 12 AS DISTANCE READING COUNTER 

L 12,-F'O 1 

* WE BEGIN TO PROCESS THE TAPE RECORD 

* 

* END OF RECORD TEST 

* 

TESEND C 5,-F'2' 

BNH EXIT 

* WE NOW ENTER THE TAPE READ LOGIC 

* 

TESTCH TM 0(3),X'F0' 

BNO DISTCE 

TM 2(3),X'F0' 

BNO ERROR 

* IF THE CHARACTER IS NUMERIC IT IS AN INTENSITY READING 

* WE CONVERT IT TO BINARY AND STORE IT IN INTTAB 

PACK DBL,0(3,3) 

* IF THE INTENSITY IS NEGITIVE WE HAVE A TAPE FORMAT ERROR 

CP DBL.-P'O' 

BNH ERROR 

CVB 7.DBL 

ST 7,0(0,*») 

* INCREMENT THE COUNTERS AND DATA BASE REGISTERS 

A 3,-F'3‘ 

A 4,=F»4' 

A 8,-F'l' 

A 6,-F'l' 

* SUBTRACT THREE FROM TOTAL BYTE COUNTER 

S 5,-F'3‘ 

B TESEND 

* IF THE CHARACTER IS ALPHABETIC IT IS A DISTANCE READING 

* WE CONVERT IT TO BINARY AND STORE IT IN D I STAB 

DISTCE TM 5(3),X'DO' 
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* 

BNO ERROR 

PACK DBL,0(6,3) 

* CONVERT THE DISTANCE READING TO A LEGITIMATE POSITIVE NUMBER 

01 DBL+7,X'0F' 

CVB 7.DBL 

ST 7,0(0,10) 

* ADVANCE THE D I STAB BASE REGISTER AND ADD 1 TO THE DISTANCE READING 

* COUNTER, SUBTRACT SIX FROM THE TOTAL BYTE COUNTER 

A 10,=F'4' 

A 12, =F' 1 1 

A 3 , =F 1 6 1 

S 5,=F'6' 

* STORE THE NUMBER OF OPTICAL DENSITY READINGS FOR THIS PEAK AND 

* ADVANCE THE PNTPER BASE REGISTER 

ST 8, 0(0, 9) 

A 9,=F'4' 

* RE INITALIZE THE INTENSITY NUMBER COUNTER 

L 8,=F'0' 

B TESEND 

* IF TAPE FORMAT IS INCORRECT GO TO ERROR EXIT 

* 


EXIT 

TM 

0(3) ,X 1 FO 


BO 

ERROR 


TM 

5 (3) ,X 1 FO 


BO 

ERROR 


* FOR CORRECT TAPE FORMAT TAKE NORMAL EXIT 

* 

* STORE FINAL DISTANCE READING IN D I STAB 

PACK DBL,0(6,3) 

01 DBL+7,X'0F' 

CVB 7,DBL 
ST 7,0(0,10) 

* INCREMENT DISTANCE READING COUNTER 

A 12, "F 1 1 1 

* STORE THE NUMBER OF OPTICAL DENSITY READINGS IN THE LAST RECORD 

* IN PNTPER 

ST 8, 0(0, 9) 

* STORE THE NUMBER OF PEAKS IN THIS RECORD IN PEKNUM 

ST 12,0(0,11) 

* WE ADD FOUR TO THE ADDRESS OF PEKNUM TO GET TO TOTP 

A 1 1 ,=F 1 4 ' 

* STORE THE TOTAL NUMBER OF OPTICAL DENSITY READINGS IN THIS RECORD 

* IN TOTP 

ST 6,0(0,11) 

* WE ADD FOUR TO THE ADDRESS OF TOTP TO GET TO THE ERROR FLAG 

A 1 1 ,=F 1 4 1 

* THIS IS A NORMAL EXIT STORE ZERO IN ERROR FLAG 

L 4,*F'0' 

ST 4,0(0,11) 

B OUT 

* WE ADD EIGHT TO THE ADDRESS OF PEKNUM TO GET TO THE ERROR FLAG 

ERROR A 1 1 ,=F ' 8 1 

* THIS IS AN ABNORMAL EXIT STORE ONE IN THE ERROR FLAG 

L 4,-F'P 

ST 4,0(0,11) 

B OUT 

* RETURN CONTROL TO CALLING PROGRAM 


OUT 

RETURN (14,12) 


BR 

14 

DBL 

DS 

D 

C0M1 

DC 

V(AMESTP) 

COM2 

DC 

V(MASPEC) 


END 

TAPE 
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C OCTOBER 10,1972 MICHAEL ROMIEZ 

C 

C THIS IS THE PROGRAM WHICH CALLS THE VARIOUS SUBROUTINES TO REDUCE THE HIGH 
C RESOLUTION MASS SPECTRAL DATA. 

C 

c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

CALL PARLST 
CALL PKLIST 
CALL STDFIT 
CALL ELELST 
STOP 
END 
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JANUARY 30,1973 

THIS SUBROUTINE IS A GENERAL POLYNOMIAL FITTING ALGORITHM. IT WILL FIT UP TO 200 
DATA POINTS TO A FROM FIRST TO TENTH DEGREE POLYNOMIAL. THE POLYNOMIAL IS BEST FIT 
IN THE LEAST SQUARES SENSE. 


X (200) 

Y (200) 

C (1 1 ) 

N 

IDEG 

SUMX (20) 

SUMXY(IO) 

SUMY 
A ( 1 1 , 1 1 ) 

B (1 1 ) 

DM ( 1 1 ) 


VARIABLES 

UP TO 200 ABCISSA VALUES 
UP TO 200 ORDINATE VALUES 

UP TO 11 COEFFICIENTS CALCULATED BY THIS SUBROUTINE 
THE NUMBER OF DATA POINTS IN THE X AND Y ARRAYS 
THE DEGREE OF THE POLYNOMIAL 

THE UP TO TWENTY VALUES OF THE SUMMATIONS OF X RAISED TO THE 
POWERS OF ONE THROUGH TWICE IDEG INCLUSIVE. 

THE UP TO TEN VALUES OF THE SUMMATIONS OF Y TIMES X RAISED 
TO THE POWERS OF ONE THROUGH IDEG INCLUSIVE. 

THE SUMMATION OF THE Y VALUES 

THE IDEG+1 BY IDEG+1 SQUARE MATRIX OF SUMMATION VALUES FROM 
THE NORMAL EQUATIONS OF THE LEAST SQUARE SYSTEM OF EQUATIONS 
THE UP TO ELEVEN CONSTANT TERMS OF THE NORMAL EQUATIONS 
TEMP STORAGE 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
SUBROUTINE POLYFT 

DOUBLE PRECISION B ( 1 1 ) ,A(1 1 , 1 1 ) ,C ( 1 1 ) ,DM ( 1 1 ) 

DOUBLE PRECISION SUMX(20) ,SUMXY(10) ,SUMY,X(200) ,Y(200) 

COMMON/FIT/X,Y,N, IDEG,C 
CLEAR VARIABLES 

DO 1000 J=1 ,10 

1000 SUMXY (j)=0. 

DO 1001 J=1 ,11 
B (J)=0. 

dm(j)=o. 

1001 C(J)=0. 

DO 1002 J=1 ,20 

1002 SUMX (J)=0. 

DO 1003 J— 1 ,11 
DO 1003 K=1 ,11 

1003 A (J ,K)=0. 

SUMY=0. 

CALCULATE THE SUMMATION VALUES 
DO 2000 M=1 , IDEG 
DO 2000 K-l.N 

2000 SUMXY (M)=SUMXY (M)+Y (K) *X (K)**M 
DO 2001 M=1 , N 

2001 SUMY=SUMY+Y(M) 

N2=l DEG*2 


DO 2002 M=1 ,N2 
DO 2002 K=1 ,N 

2002 SUMX(M)=SUMX (M)+X (K)**M 

C THIS SECTION ASSIGNS THE SUMMATION VALUES TO THE PROPER 
C MATRIX ELEMENTS 
LIMIT=IDEG+1 
DO 3001 J=l, LIMIT 
DO 3000 K=l, LIMIT 
M=J+K-2 

IF(M.EQ.O) GO TO 3000 
A(J,K)=SUMX(M) 

3000 CONTINUE 
IF(J.EQ.l) GO TO 3001 
B (J)=SUMXY (J-l ) 

3001 CONTINUE 
A ( 1 , 1 )=N 
B ( 1 )=SUMY 

C THIS SECTION CALCULATES THE ELEMENTS OF THE TRIANGULAR MATRIX 


NN=0 

DO 4001 J=1 , IDEG 
NN=NN+1 

DO 4001 K=NN , I DEG 
dm(k)=a(k+i ,j)/a(j,j> 

DO 4000 L=l, LIMIT 

4000 A(K+1 , L ) =A ( K+ 1 ,L)-DM(K)*A(J,L) 

B (K+l )=B (K+l ) -DM (K)*B (J ) 

4001 CONTINUE 

C THIS SECTION CALCULATES THE COEFFICIENTS OF THE CURVE FITTING 
C POLYNOMIAL FROM THE TRIANGULAR MATRIX 
M=LIMIT+1 
DO 5001 K=l, LIMIT 
DM(l )=0.D0 
I ROW=M-K 


I COL=M-K 

5000 IC0L=IC0L+1 



c 

c 
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IF(ICOL.GT. LIMIT) GO TO 5001 
DM( 1 )=DM ( 1 ) —A ( I ROW, I COL)*C ( I COL) 

GO TO 5000 

5001 C ( I ROW)- (B( I ROW)+DM (1 ) ) /A ( I ROW, l ROW) 
RETURN 
END 


FEBRUARY 5, 1973 


MICHAEL ROMIEZ 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


THE VARIABLES HAVE DEFALT VALUES GIVEN AT COMPILATION, BUT THEY MAY BE REVIEWED AND 
CHANGED AT EXECUTION TIME. 


VARIABLES 


DERDEG - THE DEGREE OF THE FIRST DERIVATIVE POLYNOMIAL 

DERNUM - THE NUMBER OF DATA POINTS USED BY THE FIRST DERIVATIVE POLYNOMIAL 

SMNUM - THE NUMBER OF DATA POINTS USED BY THE SMOOTHING POLYNOMIAL 

SMDEG - THE DEGREE OF THE SMOOTHING POLYNOMIAL 

CONNUM - THE NUMBER OF POINTS USED BY THE DISTANCE TO MASS CONVERSION POLYNOMIAL 

CONDEG - THE DEGREE OF THE DISTANCE OF MASS CONVERSION POLYNOMIAL 

CUTOFF - IF THE TOTAL NUMBER OF DATA POINTS IN THE PEAK PROFILE IS LESS THAN OR 

- EQUAL THIS NUMBER, THE PEAK IS TREATED AS A SINGLE PEAK WITH NO 

- DECONVOLUTION ANALYSIS 

INTMIN - IF THE PEAK PROFILE CONTAINS NO INTENSITY READINGS ABOVE THIS VALUE, THE 

- PEAK IS IGNORED. THIS ABILITY IS USEFUL IN DISCRIMINATION AGAINST SUCH 

- ARTIFACTS AS SCRATCHES ON THE PHOTOPLATE 

TOLFAC - THIS ISA FACTOR USED BY THE STANDARD MASS SEARCHING ROUTINE TO CALCULATE 

- THE TOLERANCE WITHIN WHICH IT WILL ACCEPT A PEAK AS A STANDARD MASS. 

EXPTOL - LOWER LIMIT OF THE NUMBER OF INTENSITY VALUES CONSIDERED AN 

- OVEREXPOSURE. 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
SUBROUTINE PARLST 

INTEGER DERDEG, DERNUM, SMNUM, SMDEG, CUTOFF, CONNUM, CONDEG, TOLFAC, EXPTOL 
INTEGER YES/'Y'/ 

COMMON/PAR/DERDEG , DERNUM , SMDEG , SMNUM , CUTOFF , INTMIN, CONNUM , CONDEG , TOLFAC , EXPTOL 

DERDEG=3 

DERNUM=5 

SMNUM=25 

SMDEG=3 

CONNUM=4 

C0NDEG=2 

CUT0FF=30 

I NTM I N=90 

TOLFAC=7 

EXPTOL=l 80 

WRITE(6,100) 

100 FORMAT (IX 1 DO YOU WANT TO REVIEW THE PEAK POSITION PARAMETERS?') 

READ (6,101) I 

101 FORMAT (Al) 

I F ( I .NE.YES) RETURN 
WR ITE (6, 1 02) DERDEG 

102 FORMAT (IX, 'THE DEGREE OF THE FIRST DERIVATIVE POLYNOMIAL I S', 13) 

CALL INTGER(M) 

IF(M.NE.O) DERDEG=M 
WRITE (6, 1 03) DERNUM 

103 F0RMAT(1X'THE NUMBER OF POINTS USED FOR THE FIRST DERIVATIVE IS', 13) 

CALL INTGER(M) 

IF(M.NE.O) DERNUM=M 
WRITE (6, 1 04) SMNUM 

104 FORMAT (IX, 'THE NUMBER OF POINTS USED FOR THE SMOOTHING POLYNOMIAL IS', 13) 

CALL INTGER(M) 

IF(M.NE.O) SMNUM=M 
WR ITE (6, 1 05) SMDEG 

105 FORMAT (IX, 'THE DEGREE OF THE POLYNOMIAL USED FOR SMOOTHING IS' 13) 

CALL INTGER(M) 

I f(m.ne.o) smdeg=m 
WRITE(6, 108) CONNUM 

108 FORMAT (IX, 'THE NUMBER OF MASSES USED BY THE MASS CALCULATION POLYNOMIAL IS', 13) 

CALL INTGER(M) 

IF(M.NE.O) CONNUM=M 
WRITE (6, 1 09) CONDEG 

109 FORMAT (IX,' THE DEGREE OF THE MASS CALCULATION POLYNOMIAL IS', 13) 

CALL INTGER(M) 

IF(M.NE.O) CONDEG=M 
WRITE (6, 106) CUTOFF 

106 FORMAT (IX, 'THE CUTOFF VALUE FOR A SINGLE PEAK IS',l4) 

CALL INTGER(M) 

if(m.ne.o) cutoff-m 

WRITE (6, 1 07) INTMIN 

107 FORMAT (IX, 'THE MINIMUM HIGH INTENSITY FOR A REAL PEAK IS', 14) 

CALL INTGER(M) 

IF(M.NE.O) I NTM I N=M 
: WRITE (6, 110) TOLFAC 

110 FORMAT (IX, 'THE SLIDING TOLERANCE FACTOR FOR STANDARD MASS IDENTIFICATION IS', 13) 

CALL INTGER(M) 

IF(M.NE.O) TOLFAC-M 
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C FEBRUARY 5, 1973 MICHAEL ROMIEZ 

C 

WR ITE (6, 1 1 1 ) EXPTOL 

111 FORMAT (IX, 'THE NUMBER OF INTENSITY VALUES ABOVE WHICH A PEAK WILL BE EDITED IS* , 1 4) 
CALL INTGER(M) 

IF(M.NE.O) EXPTOL=M 

RETURN 

END 
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THIS SUBROUTINE SENDS RAW PEAK PROFILES TO SUBROUTINE PEAK WHICH IN TURN CONTROLS 
THE SMOOTHING, DECONVOLUTION, AND PEAK CENTER CALCULATION ROUTINES. THIS MODULE 
SENDS UP TO 400 DATA POINTS TO SUBROUTINE PEAK AND IS RETURNED UP TO TEN PEAK CENTERS 
AND MAXIMUM INTENSITY VALUES. THESE CENTERS AND INTENSITIES ARE STORED FOR THE REST 
OF THE DATA REDUCTION ROUTINES. 


RAW I NT (25000) 
NUMINT ( 1 000 ) 

RAWDIS(IOOO) 
CENTER ( 1 500) 

H IGHT ( 1 500) 
NUMPK 
LI ST (400) 

NUM 


CENT (10) 

HGT(IO) 
LSTPAR (10) 
PTINT 
PTDIS 
PTCENT 
I NTM I N 


- TABLE OF INTENSITY VALUES FOR THE ENTIRE RAW DATA FILE 

- TABLE OF THE NUMBER OF INTENSITY READINGS PER PEAK IN THE ENTIRE 

- RAW DATA FILE 

- TABLE OF DISTANCE READINGS. ONE FOR EACH PEAK IN NUMINT 

- DISTANCE OF PEAK CENTERS FROM THE LEFT EDGE OF THE PHOTOPLATE 

- IN MILLIMETERS 

- MAXIMUM INTENSITY VALUE OF EACH PEAK 

- NUMBER OF PEAKS IN THE PEAK LIST TABLE 

- HOLDING TABLE TO PASS ONE PEAK PROFILE BETWEEN MODULES 

- EITHER THE NUMBER OF INTENSITY VALUES IN THE PEAK PROFILE OR THE 

- NUMBER OF PEAK CENTERS RETURNED BY SUBROUTINE PEAK 

- THE DISTANCE IN HALF-MICRONS OF THE LAST INTENSITY READING IN THE 

- PEAK PROFILE MEASURED FROM THE LEFT END OF THE PHOTOPLATE 

- UP TO 10 PEAK CENTERS RETURNED TO THIS SUBROUTINE FOR EACH PEAK 

- PROFILE 

- UP TO TEN MAXIMUM INTENSITY VALUES, ONE FOR EACH PEAK IN CENT 

- LIST OF PARAMETERS FROM THE PARLIST MODULE 

- POINTS TO NEXT INTENSITY READING IN RAW I NT 

- POINTS TO NEXT POSITION IN THE RAWD I S AND NUMINT ARRAYS 

- POINTS TO NEXT POSITION IN THE CENTER AND HIGHT ARRAYS 

- ALL INTENSITY READINGS IN A PROFILE MUST BE ABOVE THIS VALUE 

- OR THE PEAK IS IGNORED. THIS IS USEFUL IN DISCRIMINATING AGAINST 

- SCRATCHES ON THE PHOTOPLATE 

- LOWER LIMIT OF THE NUMBER OF INTENSITY READINGS FOR EDITING 


C XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
SUBROUTINE PKLIST 
DOUBLE PRECISION CENTER, CENT ( 1 0) 

INTEGER RAW I NT (25000) , NUM I NT (1000) ,RAWD I S ( 1000) 

INTEGER PTDIS, PTINT, LI ST (400) ,PTCENT,DIS 
INTEGER YES/'Y'/ 

DIMENSION HGT(IO) 

COMMON/PKLST/ CENTER(1500) ,HIGHT(1500) , NUMPK 
COMMON/PK/L I ST , NUM, D I S , CENT , HGT 
COMMON/PAR/LSTPAR ( 1 0) 

C SET THE I NTM IN VALUE AS GIVEN BY THE PARAMETER LIST MODULE 
INTMIN=LSTPAR(6) 

C SET THE LOWER LIMIT OF THE NUMBER OF INTENSITY READINGS FOR PROFILE EDITING 
I ED IT=LSTPAR ( 1 0) 

C CLEAR VARIABLES 

DO 100 M=1 ,1500 
CENTER(M)=0. 

100 H I GHT (M)=0. 

DO 150 M-l, 25000 

150 RAW I NT (M)=0 

DO 151 M-l, 1000 
RAWD I S (M)=0 

151 NUM I NT (M)=0 

C READ IN RAW DATA FROM DISK FILE 

READ (9) N, (RAWINT(I), 1 = 1 ,N) ,M, (NUM I NT ( I ) , RAWD I S ( I ) , 1 = 1 ,M) 

PT I NT=0 
PTD I S=1 
PTCENT=0 

10 DO 15 1-1,10 
HGT ( I )=0. 

15 CENT ( I )=0. 

C RETURN IF THERE ARE NO MORE PROFILES IN THE RAW DATA FILE 
I F (NUMINT (PTD IS). EQ. 0) GO TO 9999 

C IF THERE ARE MORE THAN 400 POINTS IN A PARTICULAR PEAK PROFILE, THE PEAK IS 
C IGNORED. A PEAK 200 MICRONS WIDE IS TOO WIDE TO BE A REAL PEAK. 

I F (NUM I NT (PTD IS). GT. 400) GOTO 11 
I F (NUM I NT (PTD I S) .GT. 4) GO TO 200 

C IF A PROFILE HAS FEWER THAN 5 INTENSITY READINGS IT IS DELETED AS A NOISE SPIKE 

11 PTD I S=PTD I S+l 

PT I NT=PT I NT+NUM I NT (PTD I S- 1 ) 

GO TO 10 

C READ PEAK PROFILE AND DISTANCE INTO DATA PASSING ARRAYS 
200 NUM=NUM I NT ( PTD I S ) 

DO 225 1-1,400 
225 LI ST( I )=0 

DO 250 J=1 ,NUM 
PT I NT=PT I NT+1 
250 L I ST (J) -RAW I NT (PTINT) 

D I S=RAWD I S ( PTD I S ) 
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C JANUARY 11 ,1973 
C 

C IF THE HIGHEST INTENSITY VALUE IS LESS THAN INTMIN IGNORE THE PEAK 
DO 275 J-l.NUM 

IF(LIST(J).GT. INTMIN) GO TO 290 
275 CONTINUE 

PTDIS-PTDIS+1 
GO TO 10 

C IF THERE ARE MORE THAN I ED IT INTENSITY READINGS IN A PEAK PROFILE, IT IS AN 
C OVEREXPOSURE AND IS DISPLAYED FOR POSSIBLE EDITING. 

290 I F (NUM.GT. (EDIT) CALL EDIT 
C SEND PEAK PROFILE TO BE ANALYZED 
CALL PEAK 

C KEEP PEAK CENTERS AND INTENSITIES IN CENTER AND HIGHT ARRAYS 
DO 299 M-1.NUM 
PTCENT-PTCENT+1 

C THE PEAK CENTER VALUES ARE CONVERTED FROM HALF-MICRONS TO MILLIMETERS 
CENTER(PTCENT)«CENT(M)/2000. 

C THE DISTANCE ENCODER RECYCLES AT 255 MILLIMETERS AND THIS TEST 
C COMPENSATES FOR RECYCLING. 

IF(PTCENT.EQ. 1) GO TO 299 

I F (CENTER (PTCENT) .LT. CENTER (PTC ENT- I ) .AND .CENTER (PTC ENT- 1 ) - 
X.GT.200) CENTER (PTC ENT) -CENTER (PTC ENT) +255- 
299 HE I GHT (PTCENT)“HGT (M) 

C THIS TEST WILL TRUNCATE THE PEAK LIST IF IT APPROACHES WITHIN 10 VARIABLES OF 
C THE END OF THE ARRAY 

I F (PTCENT. GT. 1490) GO TO 400 
PTDIS-PTDIS+1 
GO TO 10 
400 WRITE(6,1) 

1 FORMAT ( 'O' , 'THE PEAK LIST WAS TRUNCATED AT', 1 6,' VALUES') 

9999 NUMPK-PTCENT 
RETURN 
END 


MICHAEL ROM I EZ 


C JANUARY 17,1973 
C 

C THIS SUBROUTINE CONTROLS DATA SMOOTHING, PEAK DECONVOLUTION, AND PEAK CENTER CALCULATION 
C 

C VARIABLES 

C 

C LI ST (400) - HOLDING TABLE TO PASS ONE PEAK PROFILE BETWEEN MODULES 

C NUM - EITHER THE NUMBER OF INTENSITY VALUES IN THE PEAK PROFILE OR THE 

C - NUMBER OF PEAK CENTERS RETURNED BY SUBROUTINE ANAL 

C DIS - THE DISTANCE IN HALF-MICRONS OF THE LAST INTENSITY READING IN THE 

C - PEAK PROFILE MEASURED FROM THE LEFT END OF THE PHOTOPLATE 

C CENT (10) - UP TO 10 PEAK CENTERS RETURNED TO THIS SUBROUTINE FOR EACH PEAK 

C - PROFILE 

C HGT(IO) - UP TO TEN MAXIMUM INTENSITY VALUES, ONE FOR EACH PEAK IN CENT 

C LSTPAR(IO) - LIST OF PARAMETERS FROM THE PARLIST MODULE 

C NUMSM - NUMBER OF POINTS USED BY THE SMOOTHING ROUTINE 

C CUTOFF - MINIMUM NUMBER OF POINTS FOR A PROFILE TO BE ANALYZED FOR 

C - MORE THAN ONE PEAK. PROFILES WITH LESS THAN CUTOFF DATA 

C - POINTS WILL BE TREATED AS A SINGLE PEAK 

C 

c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE PEAK 

DOUBLE PRECISION CENT ( 1 0) 

INTEGER DIS, CUTOFF 
DIMENSION HGT (10) .LIST (400) 

COMMON/PAR/LSTPAR ( 1 0) 

COMMON/PK/L I ST , NUM , D I S , CENT , HGT 

C SET THE PARAMETERS FOR THE NUMBER OF POINTS USED BY THE SMOOTHING ROUTINE AND 
C THE MINIMUM NUMBER OF POINTS IN A PROFILE CONSIDERED FOR DECONVOLUTION ANALYSIS 
NUMSM=LSTPAR(A) 

CUT0FF=LSTPAR(5) 

C IF THERE ARE LESS THAN CUTOFF DATA POINTS, THIS IS A SINGLE PEAK 
IF(NUM.GT. CUTOFF) GO TO 100 

C IF THERE ARE ENOUGH RAW DATA POINTS IN THE PROFILE, SMOOTH THE DATA 
C BEFORE SENDING IT TO MODULE ONEPK 

I F (NUM.GT. (NUMSM+1 0) ) GO TO 150 

CALL ONEPK 

RETURN 

150 CALL SMOOTH 

C THE DATA SMOOTHING ROUTINE RETURNS (NUMSM- 1) LESS DATA POINTS THAN WAS IN THE 
C RAW DATA. THE DISTANCE VALUE WHICH INDICATES THE DISTANCE OF THE LAST DATA POINT 
C IN THE LIST MUST BE ADJUSTED 
D I S=D I S- ( (NUMSM- 1 ) /2) 

CALL ONEPK 
RETURN 

C SMOOTH THE DATA 
100 CALL SMOOTH 

C THE DATA SMOOTHING ROUTINE RETURNS (NUMSM- 1) LESS DATA POINTS THAN WAS IN THE 
C RAW DATA. THE DISTANCE VALUE WHICH INDICATES THE DISTANCE OF THE LAST DATA POINT 
C IN THE LIST MUST BE ADJUSTED 
D I S=D I S- ( (NUMSM- 1 ) /2) 

C SUBROUTINE ANAL USES THE FIRST DERIVATIVE TO DETERMINE IF THERE IS MORE THAN ONE 
C PEAK IN A PARTICULAR SMOOTHED PROFILE 
CALL ANAL 
RETURN 
END 
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C JANUARY 16,1973 

C 

C THIS ROUTINE CALLS A GENERAL LEAST SQUARE POLYNOMIAL FITTING ROUTINE TO SMOOTH 

C THE PROFILE DATA. IT RETURNS NUM LESS (NUMSUM-1) SMOOTHED DATA POINTS. THE DEGREE 

C OF THE POLYNOMIAL AND THE NUMBER OF DATA POINTS USED FOR SMOOTHING IS SELECTED IN 

C PARAMETER SELECTING SUBROUTINE PARLST. 

C 

C VARIABLES 

C 

C L I ST (400) - HOLDING TABLE TO PASS ONE PEAK PROFILE BETWEEN MODULES 

C NUM - THE NUMBER OF INTENSITY VALUES IN THE PEAK PROFILE 

C DIS - THE DISTANCE IN HALF-MICRONS OF THE LAST INTENSITY READING IN THE 

C - PEAK PROFILE MEASURED FROM THE LEFT END OF THE PHOTOPLATE 

C CENT (10) - UP TO 10 PEAK CENTERS 

C HGT(IO) - UP TO TEN MAXIMUM INTENSITY VALUES, ONE FOR EACH PEAK IN CENT 

C LSTPAR (10) - LIST OF PARAMETERS FROM THE PARLIST MODULE 

C X (200) - UP TO 200 ABCISSA VALUES TO BE PASSED TO THE POLYNOMIAL FITTING 

C - ROUTINE 

C Y (200) - UP TO 200 ORDINATE VALUES TO BE PASSED TO THE POLYNOMIAL FITTING 

C - ROUTINE 

C N - NUMBER OF X AND Y VALUES PASSED 

C I DEG - DEGREE OF THE POLYNOMIAL (UP TO TENTH POSSIBLE) 

C C (1 1 ) - UP TO 11 COEFFICIENTS RETURNED BY THE POLYNOMIAL ROUTINE 

C DEGSM - DEGREE OF THE SMOOTHING POLYNOMIAL, SUPPLIED BY THE PARLST MODULE 

C NUMSM - NUMBER OF DATA POINTS USED BY THE SMOOTHING POLYNOMIAL, SUPPLIED 

C - BY THE PARLST MODULE 

C HOLD (400) - TEMPORARY STORAGE USED BY THIS MODULE 

C 

C XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
SUBROUTINE SMOOTH 

DOUBLE PRECISION X(200) ,Y(200) ,C (1 1) ,CENT(10) 

DIMENSION LIST(400) ,HGT (10) 

INTEGER DIS, DEGSM, H0LD(400) 

COMMON/PAR/LSTPAR ( 1 0) 

COMMON/F I T/X , Y , N , I DEG , C 
COMMON/PK/L I ST , NUM , D I S , CENT , HGT 

C SET THE DEGREE AND NUMBER OF POINTS USED BY THE SMOOTHING POLYNOMIAL 
DEGSM=LSTPAR(3) 

NUMSM=LSTPAR(4) 

C CLEAR TEMP STORAGE 
DO 100 M=1 ,400 
100 HOLD (M)=0 

L IM !T=NUMSM- 1 
M-NUM-L IM IT 

C LOAD RAW DATA POINTS INTO X AND Y ARRAYS 
DO 10 I *2, NUMSM 
J= I - 1 

10 Y ( I J-DFLOAT ( L I ST ( J ) ) 

DO 150 L-l, NUMSM 

150 X (L)=DFLOAT (L) 

I DEG-DEGSM 
N-NUMSM 

C THIS DO LOOP (DOWN TO STMT 200) CONTROLS THE SENDING OF NUMSM RAW DATA POINTS TO THE 
C POLYNOMIAL ROUTINE TO OBTAIN THE SMOOTHED DATA POINTS. ALL NUMSM DATA POINTS ARE USED 
C TO CALCULATE ONE SMOOTH DATA POINT. THE LOWEST RAW DATA POINT IS DISCARDED AND THE 
C REMAINING POINTS ARE MOVED DOWN ONE POSITION. A NEW HIGH DATA POINT IS READ IN FROM THE 
C PROFILE DATA AND THE NEXT SMOOTHED DATA POINT IS CALCULATED. THE PROCESS IS REPEATED 
C UNTIL WE HAVE NUM- (NUMSM- 1 ) SMOOTHED DATA POINTS. 

DO 200 1-1 ,M 
J-I+LIMIT 
DO 11 K-l, LIMIT 
KA-K+1 

11 Y(K)-Y(KA) 

Y (NUMSM) -DFLOAT ( L I ST ( J ) ) 

CALL POLYFT 
I NDEX1-J-L IM IT 
I NDEX2-LIMIT/2+1 
HOLD ( I NDEX1 )-C ( 1 ) 

DO 200 L-l, DEGSM 
LL-L+1 

C PUT THE CALCULATED SMOOTH DATA POINT INTO ARRAY HOLD FOR TEMP STORAGE 
200 HOLD ( INDEX1 )-HOLD ( I NDEX 1 )+C (LL)*FLOAT ( I NDEX2 ) **L 
C CHANGE THE VALUE OF NUM TO THE NEW NUMBER OF DATA POINTS IN THE PROFILE AFTER 
C SMOOTHING 

NUM-NUM-L IM IT 

C THE LAST TWO DO LOOPS MOVE THE SMOOTH DATA POINTS FROM TEMP STORAGE BACK TO THE 
C PROFILE PASSING ARRAY 
DO 300 L-l, NUM 

300 LI ST (L)-HOLD (L) 

I STRT-NUM+1 

I STOP-NUM+NUMSM 
DO 301 K-ISTRT, I STOP 

301 L I ST (K)»0 
RETURN 
END 
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JANUARY 16,1973 

THIS MODULE CALLS THE GENERAL LEAST SQUARE FITTING ROUTINE TO CALCULATE A FIRST DERIVATIVE 
OF THE SMOOTHED PEAK PROFILE. THE FIRST DERIVATIVE IS FOUND BY OVERLAPPING POLYNOMIALS 
IN THE SAME MANNER AS THE SMOOTHING ALGORITHM. 


LIST(400) 

NUM 

DIS 

CENT (10) 

HGT(IO) 

LSTPAR(IO) 

X (200) 

Y (200) 

N 

I DEG 
C(ll) 

DEGFT 

NUMFT 

FSTDER(400) 

INDEX1 

INDEX2 

INDEX3 


VARIABLES 

HOLDING TABLE TO PASS ONE PEAK PROFILE BETWEEN MODULES 

THE NUMBER OF INTENSITY VALUES IN THE PEAK PROFILE OR THE NUMBER 

OF PEAKS FOUND BY THIS MODULE 

THE DISTANCE IN HALF-MICRONS OF THE LAST INTENSITY READING IN THE 
PEAK PROFILE MEASURED FROM THE LEFT END OF THE PHOTOPLATE 
UP TO 10 PEAK CENTERS 

UP TO TEN MAXIMUM INTENSITY VALUES, ONE FOR EACH PEAK IN CENT 
LIST OF PARAMETERS FROM THE PARLIST MODULE 

UP TO 200 ABCISSA VALUES TO BE PASSED TO THE POLYNOMIAL FITTING 
ROUTINE 

UP TO 200 ORDINATE VALUES TO BE PASSED TO THE POLYNOMIAL FITTING 
ROUTINE 

NUMBER OF X AND Y VALUES PASSED 

DEGREE OF THE POLYNOMIAL (UP TO TENTH POSSIBLE) 

UP TO 11 COEFFICIENTS RETURNED BY THE POLYNOMIAL ROUTINE 
DEGREE OF THE 1ST DERIVATIVE POLYNOMIAL, SUPPLI ED BY THE PARLST MODULE 
NUMBER OF DATA POINTS USED BY THE 1ST DERIVATIVE POLYNOMIAL, SUPPLIED 
BY THE PARLST MODULE 

TABLE OF FIRST DERIVATIVE POINTS AS CALCULATED BY THIS MODULE 
POINTER FOR THE FSTDER ARRAY 

POINTS TO THE CENTER OF THE NUMFT POINTS USED BY THE POLYNOMIAL 
POINTS TO THE LAST USED POSITION IN THE FSTDER ARRAY 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
SUBROUTINE ANAL 

DOUBLE PRECISION CENT(IO) ,C(1 1 ) ,X(200) ,Y(200) 

DIMENSION L I ST(400) ,HGT(10) ,FSTDER(400) 

INTEGER DIS, DEGFT, HOLD (400) 

COMMON/F IT/X, Y,N, IDEG,C 
COMMON/PAR/LSTPAR (10) 

COMMON/PK/L I ST , NUM , D I S , CENT , HGT 

SET THE DEGREE AND NUMBER OF POINTS USED BY THE 1ST DERIVATIVE POLYNOMIAL 
DEGFT=LSTPAR(1 ) 

NUMFT=LSTPAR(2) 

CLEAR VARIABLES 

DO 100 M=1 ,400 
100 FSTDER(M)=0. 

DO 110 M=1 ,200 

X (m)=o. 

110 Y (M)=0. 


I DEG=DEGFT 
N=NUMFT 
LIMIT-NUMFT-1 
M=NUM-L IM IT 

C LOAD DATA POINTS INTO X AND Y ARRAYS 
DO 10 I -2, NUMFT 
J= I - 1 

10 Y ( I )=DFLOAT (L I ST (J) ) 

DO 150 L-l f NUMFT 

150 X (L)=DFLOAT (L) 

C THIS DO LOOP (DOWN TO STMT 200) CONTROLS THE SENDING OF NUMFT DATA POINTS TO THE 
C POLYNOMIAL ROUTINE TO OBTAIN THE 1ST DERIVATIVE CURVE. ALL NUMSM DATA POINTS ARE USED 
C TO CALCULATE ONE 1ST DERIVATIVE DATA POINT. THE LOWEST DATA POINT IS DISCARDED AND THE 
C REMAINING POINTS ARE MOVED DOWN ONE POSITION. A NEW HIGH DATA POINT IS READ IN FROM THE 
C PROFILE DATA AND THE NEXT 1ST DERIVATIVE POINT IS CALCULATED. THE PROCESS IS REPEATED 
C UNTIL WE HAVE NUM- (NUMFT-1 ) 1ST DERIVATIVE POINTS. 

DO 200 1-1, M 
J=I+LIMIT 
DO 11 K=l, LIMIT 
KA=K+1 

11 Y (K)=Y (KA) 

Y ( NUMFT )=D FLOAT (LIST (J) ) 

CALL POLYFT 
I NDEX1 -J-L IM IT 
I NDEX2=LIM IT/2+1 
FSTDER ( I NDEX1 )=C (2) 

DO 200 L=2, DEGFT 
LL=L+1 


LLL=L- 1 

200 FSTDER ( INDEX1 )=FSTDER ( I NDEX1 )+D FLOAT (L)*C (LL)*DEFLOAT ( I NDEX2)**LLL 
INDEX3-NUM-LIMIT 

C THIS SECTION SEARCHES THROUGH THE FIRST DERIVATIVE CURVE TO LOCATE THE CENTER OF AND 
C HEIGHT OF UP TO TEN PEAKS CONTAINED IN EACH PEAK PROFILE. IN THIS SECTION J IS THE CENT 
C AND HGT ARRAY POINTER AND M IS THE FSTDER ARRAY POINTER. 



C JANUARY 16,1973 MICHAEL ROMIEZ 

C 

J=1 

M=0 

1 M-M+l 

I F (M . GT . INDEX3-1 ) GO TO 99 
C UPSLOPE LOGIC 

IF(FSTDER(M+1)-FSTDER(M))3,1 ,1 

2 M=M+ 1 

IF (M.GT. INDEX3) GO TO 99 

C WE HAVE JUST STARTED DOWN THE SIDE OF THE FIRST DERIVATIVE CURVE. WE CONTINUE DOWN 
C UNTIL WE GO THROUGH THE X-AXIS. THE POINT AT WHICH THE FIRST DERIVATIVE CURVE PASSES 
C THROUGH THE X-AXIS GOING FROM POSITIVE TO NEGATIVE, IS THE CENTER OF A PEAK. 

3 I F (FSTDER(M) .GT.O. ) GO TO 2 

C THE NEXT SIX STATEMENTS CALCULATE THE CENTER AND HEIGHT OF THE JUST FOUND PEAK 
DHIGH-FSTDER(M-l) 

DLOW-ABS (FSTDER (M) ) 

CENT ( J ) -FLOAT (M- 1 )+DH I GH/DH l GH+DLOW) 

CENT (J)-CENT (j)+LIMIT/2 

K-M+LIMIT/2 

HGT (J)-L I ST (K) 

C IF THIS MODULE HAS ALREADY FOUND TEN PEAKS, STOP LOOKING 
J-J+l 

IF(J.GT. 10) GO TO 99 

C THIS LOOP TESTS IF THE FIRST DERIVATIVE CURVE GOES BACK ABOVE THE X-AXIS. IF IT 
C DOES IT CAUSES A BRANCH TO THE UPSLOPE LOGIC SECTION 
DO 80 K-M, I NDEX3 
I F (FSTDER ( K) .LT.O.) GO TO 80 
GO TO 101 
80 CONTINUE 
GO TO 99 
101 M-K-l 
GO TO 1 
99 NN-J-1 

C THE CENTER OF EACH PEAK FOUND BY THIS MODULE IS CALCULATED IN TERMS OF THE DISTANCE 
C IN HALF-MICRONS FROM THE LEFT EDGE OF THE PHOTO PLATE. 

DO 89 K-l.NN 

89 CENT (K)-D I S- (NUM-CENT (K) ) 

C NUM IS SET TO THE NUMBER OF PEAKS FOUND 
NUM-NN 
RETURN 
END 
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THIS SUBROUTINE WILL DETERMINE THE CENTER AND HIGHT OF A SINGLE PEAK. IT FITS 
THE PEAK PROFILE TO A GAUSSIAN DISTRIBUTION CURVE. 


LIST(AOO) 

NUM 


CENT (10) 
HGT(IO) 


- UP TO 200 POINTS THAT ARE USED BY THE POLYNOMIAL FITTING 

- ROUTINE FOR THE BEST LEAST SQUARES GAUSSIAN FIT 

- ONE ORDINATE VALUE FOR EACH X VALUE 

- NUMBER OF X AND Y VALUES 

- IN THIS CASE IDEG=2. IT IS THE DEGREE OF THE POLYNOMIAL USED 

- FOR FITTING 

- THE FIRST THREE VALUES OF THE ARRAY CONTAIN THE COEFFICIENTS OF 

- THE SECOND DEGREE POLYNOMIAL 

- HOLDING TABLE TO PASS PEAK PROFILES BETWEEN MODULES 

- EITHER THE NUMBER OF INTENSITIES IN THE PEAK PROFILE OR NUM-I 

- TO INDICATE THE NUMBER OF PEAKS FOUND BY THIS MODULE 

- THE DISTANCE IN HALF-MICRONS OF THE LAST INTENSITY READING IN 

- THE PROFILE AS MEASURED FROM THE LEFT END OF THE PHOTOGRAPHIC 

- PLATE 

- THE CENTER CALCULATED BY THIS MODULE IS PUT IN CENT (1 ) 

- THE HEIGHT CALCULATED BY THIS MODULE IS PUT IN HGT(l) 


c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE ONEPK 

DOUBLE PRECISION CENT (1 0) ,C(1 1 ) ,X(200) ,Y(200) 

INTEGER DIS 

DIMENSION HGT (10) , L I ST (400) 

COMMON/FIT/X,Y,N, IDEG.C 

COMMON/PK/L I ST , NUM , D I S , CENT , HGT 

I F (NUM. GT.200) RETURN 

N-NUM 

IDEG-2 

C PUT THE LOG OF THE INTENSITY VALUE IN THE ORDINATE LIST 
DO 50 J-l.NUM 
X(J)-DFL0AT(J) 

50 Y (J)-DLOG (DFLOAT (L I ST (J) ) ) 

C CALL THE POLYNOMIAL FITTING ROUTINE 
CALL POLYFT 

I F (C (3) . GT.O. ) GO TO 99 

C CALCULATE THE CENTER AND HEIGHT OF THE PEAK FROM THE COEFFICIENTS OF THE POLYNOMIAL 
CENT (1 )*(— C (2)/C(3) )/2. 

HGT(1)-C(1) + (CENT(1)/SQRT(-1./C(3)))**2 
HGT(1)-EXP(HGT(1)) 

CENT ( 1 )-D I S- (NUM-CENT (1 ) ) 

NUM-1 
99 RETURN 
END 
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C FEBRUARY 3,1973 MICHAEL ROMIEZ 

C 

C THIS MODULE WILL ALLOW EDITING OF THE RAW PEAK PROFILE. PART OF THE PROFILE MAY BE 

C DELETED AND SELECTED INTENSITY READINGS MAY BE CHANGED. 

C 

C VARIABLES 

C 

C LIST (400) - HOLDING TABLE TO PASS ONE PEAK PROFILE BETWEEN MODULES 

C NUM - THE NUMBER OF INTENSITY VALUES IN THE PEAK PROFILE 

C DIS - THE DISTANCE IN HALF-MICRONS OF THE LAST INTENSITY READING IN THE 

C - PEAK PROFILE MEASURED FROM THE LEFT END OF THE PHOTOPLATE 

C CENT (10) - UP TO 10 PEAK CENTERS 

C HGT(IO) - UP TO TEN MAXIMUM INTENSITY VALUES, ONE FOR EACH PEAK IN CENT 

C 

c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE EDIT 
DOUBLE PRECISION CENT (10) 

INTEGER LIST (400) ,D I S 
DIMENSION HGT(IO) 

COMMON/PK/L I ST , NUM , D I S , C ENT , HGT 
DATA lYES/'Y 1 / 

C THIS SECTION LISTS THE RAW PEAK PROFILE 
WRITE (6, 1 00) 

100 F0RMAT('0' , 'PEAK PROFILE') 

WRITE(6,101)(LIST(I),I=1 ,NUM) 

101 FORMAT (5X, 20 1 5) 

C THIS SECTION WILL DELETE PART OF THE PEAK PROFILE. IT WILL ACCEPT A LOW AND HIGH 
C SUBSCRIPT AND ADJUST THE INTENSITY READINGS IN ARRAY LIST TO THESE NEW VALUES. 

C IF THE NEW LOW SUBSCRIPT ENTERED IS GREATER THAN NUM, THIS SECTION WILL BE SKIPPED. 

WRITE (6, 1 02) 

102 FORMAT ( 'O' , ' ENTER A NEW STARTING INTENSITY') 

CALL INTGER(LOW) 

IF(LOW.GT.NUM) GO TO 400 
IF (I.EO.l) GO TO 200 
L IMIT“NUM- (LOW- 1 ) 

DO 150 J-l, LIMIT 
INDEX=L0W-1+J 
150 L I ST(J)=LI ST ( I NDEX) 

200 WR ITE (6, 1 04) 

104 FORMAT (' O' , 'ENTER NEW ENDING SUBSCRIPT') 

C IF NOTHING IS ENTERED FOR THE ENDING SUBSCRIPT, THE OLD VALUE IS KEPT. 

CALL I NTGER( I HIGH) 

IF(IHIGH.EQ.O) IHIGH=NUM 
IF(IHIGH.EQ.NUM) GO TO 300 
ISTRT=IHIGH+1 
DO 250 J=l STRT.NUM 
250 LIST(J)=0 

D I ST=D I S- (NUM- 1 HI GH) 

300 NUM=NUM- (NUM- IHIGH) - (LOW-1 ) 

C THIS SECTION WILL ACCEPT INTENSITY VALUE CHANGES 
400 WRITE (6, 105) 

105 FORMAT ('O', 'ANY INTENSITY EDITING?') 

READ (5, 106) I 

106 FORMAT (Al) 

IF(I .NE. IYES) RETURN 
WRITE(6,101) (LIST(|),I«1 ,NUM) 

405 WRITE (6, 107) 

107 FORMAT(IX, 'ENTER SUBSCRIPT') 

CALL INTGER(M) 

IF(M.EQ.O) RETURN 
WRITE (6, 1 08) 

108 FORMAT (IX, 1 ENTER INTENSITY') 

CALL INTGER(N) 

LIST(M)=N 
GO TO 450 
END 
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JANUARY 26,1973 MICHAEL ROMIEZ 

THIS SUBROUTINE CONTROLS THE STANDARD MASS SEARCH AND LISTING ROUTINES 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


C 


C 


c 


c 

c 


c 

c 


c 

c 

c 


c 

c 


VARIABLES 


CENTER(1500) 

H I GHT ( 1 500) 
NUMPK 

STDD I S (160) 
STDMAS (160) 
NUMSTD 
REFMAS ( 1 60) 


- DISTANCE OF PEAK CENTERS FROM THE LEFT EDGE OF THE PHOTOPLATE 

- IN MILLIMETERS 

- MAXIMUM INTENSITY VALUE OF EACH PEAK 

- NUMBER OF PEAKS IN THE PEAK LIST TABLE 

- THE DISTANCE VALUES OF PEAK CENTERS IDENTIFIED AS STANDARD MASSES 

- THE SQUARE ROOT OF THE MASS VALUE OF IDENTIFIED STANDARD MASSES 

- THE NUMBER OF IDENTIFIED STANDARD MASSES 

- A TABLE OF THE MOST INTENSE LINES IN THE MARKER COMPOUND SPECTRUM. 

- THESE ARE THE MASSES SEARCHED FOR BY THE STANDARD MASS IDENTIFICATION 

- ROUTINE 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE STDFIT 

DOUBLE PRECISION C ENTER, STDD I S(l60) , STDMAS (160) , REFMAS (160) 

INTEGER YES/'Y'/ 

COMMON/PKLST/ CENTER (1 500) ,H I GHT (1 500) , NUMPK 
COMMON/PFKLST/ STDD I S , STDMAS , REFMAS , NUMSTD 
CLEAR VARIABLES 

DO 99 M-l ,160 
STDD I S(M)*0. 

99 STDMAS (M)-0. 

READ IN THE STANDARD MASSES AND GET THEIR SQUARE ROOT 
READ (7) REFMAS 
DO 100 M=1 ,160 

100 REFMAS (M)=DSQRT (REFMAS (M)) 

LIST UP TO THE FIRST FORTY DISTANCE VALUES IN ARRAY CENTER ON THE OUTPUT TERMINAL 
I =NUMPK 

IF(I.GT.40) 1=40 
WR ITE (6, 1 ) I 

1 FORMAT (' O' , 'THESE ARE THE FIRST' , 13, 'DISTANCE VALUES') 

WRITE(6,2) (CENTER (J) ,J=1 ,1) 

2 FORMAT ( IX , 1 0F1 2. 5) 

NUMSTD-0 

3 NUMSTD-NUMSTD+1 

THE PROGRAM USER MUST ENTER A KNOWN MASS VALUE FOR AT LEAST TWO OF THE MASS VALUES 
IN ARRAY CENTER. 

WRITE (6,4) 

4 FORMAT (IX,' ENTER SUBSCRIPT') 

CALL INTGER(I) 

I F ( I . EQ.O) GO TO 66 
WRITE(6,5) 

5 FORMAT (IX, 'ENTER MASS') 

READ (5,6) STDMAS (NUMSTD) 

6 FORMAT (FI 0.0) 

STDD I S (NUMSTD)=CENTER ( I ) 

STDMAS (NUMSTD ) =DSQRT (STDMAS (NUMSTD ) ) 

GO TO 3 

66 NUMSTD=NUMSTD-1 

AT THIS POINT THE USER CAN TRANSFER CONTROL TO THE STANDARD MASS MANUAL ADDITION MODULE 
SO THAT HE CAN USE THE LIST EDITING CAPABILITIES OF THAT MODULE. 

WRITE(6,7) 

7 FORMATCO' , 'ANY ADDITIONS?') 

READ (5,8) I 

8 FORMAT (A 1) 

I F ( I .NE.YES) GO TO 200 
GO TO 600 

CALL THE STANDARD MASS FINDER SUBROUTINE 
200 CALL STDHNT 

LIST THE CURRENT STANDARD MASSES IDENTIFIED BY THE LAST PASS OF THE FINDER ROUTINE 
300 CALL STDLST 

AT THIS POINT THE USER MAY EDIT THE STANDARD MASS LIST 
400 WRITE(6,9) 

9 F0RMAT('0' , 'ANY CHANGES?') 

READ(5,8) I 

KEY* I 

I F ( I . EQ. YES) GO TO 500 
450 WRITE(6,10) 

10 FORMATCO' , 'ANY ADDITIONS?') 

READ(5,8) I 

I F ( I . EQ. YES) GO TO 600 

IF THERE ARE NO CHANGES OR ADDITIONS, THE CURRENT LIST OF STANDARD MASSES IS USED 
BY THE DISTANCE TO MASS CONVERSION POLYNOMIAL FOR LEAST SQUARE POLYNOMIAL FITTING. 

IF(KEY. NE.YES) RETURN 
GO TO 200 
500 CALL CHANGE 
GO TO 450 
600 CALL ADD 
GO TO 200 
END 
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THIS MODULE WILL SEARCH FOR THE MAKER COMPOUND PEAKS IN ARRAY CENTER 


CENTER(1500) 

H I GHT ( 1 500) 
NUMPK 

STDD I S (160) 
STDMAS (160) 
NUMSTD 
REFMAS(160) 


LSTPAR(IO) 

X (200) 

Y (200) 

N 

I DEG 
C(ll) 

ITOL 


K 

M 

MSTRT 

MSTOP 

KSTOP 

TEMP 


VARIABLES 

DISTANCE OF PEAK CENTERS FROM THE LEFT EDGE OF THE PHOTOPLATE 
IN MILLIMETERS 

MAXIMUM INTENSITY VALUE OF EACH PEAK 
NUMBER OF PEAKS IN THE PEAK LIST TABLE 

THE DISTANCE VALUES OF PEAK CENTERS IDENTIFIED AS STANDARD MASSES 
THE SQUARE ROOT OF THE MASS VALUE OF IDENTIFIED STANDARD MASSES 
THE NUMBER OF IDENTIFIED STANDARD MASSES 

A TABLE OF THE MOST INTENSE LINES IN THE MARKER COMPOUND SPECTRUM. 

THESE ARE THE MASSES SEARCHED FOR BY THE STANDARD MASS IDENTIFICATION 
ROUT I NE 

PARAMETER LIST FROM PARAMETER LIST MODULE 

UP TO 200 POINTS THAT ARE USED FOR POLYNOMIAL FITTING 

ONE ORDINATE VALUE FOR EACH X VALUE 

NUMBER OF X AND Y VALUES 

IT IS THE DEGREE OF THE POLYNOMIAL USED FOR FITTING 

UP TO ELEVEN COEFFICIENTS OF THE POLYNOMIAL CALCULATED 

THIS IS A FACTOR USED BY THE STANDARD MASS SEARCHING ROUTINE TO CALCULATE 

THE TOLERANCE WITHIN WHICH IT WILL ACCEPT A PEAK AS A STANDARD MASS. A SMALL 

VALUE MEANS A SMALL TOLERANCE 

THE REFMAS ARRAY POINTER 

THE CENTER ARRAY POINTER 

LOWER LIMIT FOR VALUES OF M 

UPPER LIMIT FOR VALUES OF M 

UPPER LIMIT FOR VALUES OF K 

THE DISTANCE WHERE THE NEXT STANDARD MASS IS EXPECTED 
TO BE FOUND 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
SUBROUTINE STDHNT 

DOUBLE PRECISION CENTER, STDD I S(l60) ,STDMAS(160) ,REFMAS(I60) 

DOUBLE PRECISION X (200) , Y (200) ,C(1 1 ) ,TEMP,HTOL, LTOL 
COMMON/PKLST/ CENTER ( 1 500) , H I GHT ( 1 500) , NUMPK 
COMMON/PFKLST/ STDD I S , STDMAS , REFMAS , NUMSTD 
COMMON/FIT/ X,Y,NUM, IDEG.C 
COMMON/PAR/LSTPAR ( 1 0) 

SET THE SLIDING TOLERANCE FACTOR 
IT0L»LSTPAR(9) 

CLEAR VARIABLES 

DO 100 M-l ,200 
X(M)«0. 

100 Y(M)«0. 

WE WANT TO USE A QUADRATIC SEARCH ROUTINE FOR WHICH WE NEED 
THREE PEAKS. THE FIRST TIME WE ENTER THIS SUBROUTINE WE MAY 
ONLY HAVE TWO IDENTIFIED STANDARD MASSES. A LINAER SEARCH SECTION 
IS PROVIDED FOR THIS CASE 

IF (NUMSTD. GT. 2) GO TO 200 
DO 150 M-l, 2 

x ( m)-stdd i s (m) 

Y (M) -STDMAS (M) 

NUM-2 
I DEG— 1 

CALL POLYFT 
GO TO 275 

THIS IS THE QUADRATIC SEARCH SECTION 
200 DO 250 M-l, 3 
x ( m)-stdd 1 s ( m) 

Y(M) -STDMAS (M) 

NUM-3 
IDEG-2 
CALL POLYFT 

IS THE SEARCHING SECTION 
K-0 
M-l 

KSTOP- 160 
MSTOP-NUMPK 


150 


250 


THIS 

275 


301 K-K+l 

C IF WE HAVE LOOKED FOR ALL THE REFMAS ENTRIES, RETURN 
I F (K.GT . KSTOP) RETURN 
IF(NUMSTD.LT.3) GO TO 350 

C THIS IS THE SECTION WHERE THE EXPECTED DISTANCE VALUE FOR THE NEXT ENTERY 
C IN THE REFMAS ARRAY IS CALCULATED. 

TEMP-(-C(2)+DSQRT(C(2)**2-i»*C (3)* (C ( 1 ) -REFMAS (K))))/(2*C (3)) 

GO TO 351 

350 TEMP- (REFMAS (K) -C ( 1 ) )/C (2) 

C WE MUST SPECIFY A TOLERANCE WITHIN WHICH WE WILL ACCEPT A 
C DISTANCE VALUE AS THE SEARCHED FOR STANDARD MASS. WE HAVE A 
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C 

C FORMULA WHICH ALLOWS A HIGHER TOLERANCE AT THE HIGH MASS END. IT 
C USES A FACTOR (ITOL) WHICH IS SUPPLIED BY THE PARAMETER 
C SETTING MODULE. 

351 TOL= (FLOAT ( ITOL)*. 001 / (TEMP/350. ) 

HTOL=TEMP+TOL 
LTOL=TEMP-TOL 
GO TO A02 
1+01 M=M+1 

C SEARCH ONLY WITHIN THE SUBLIST SPECIFIED BY MSTRT AND MSTOP. 

1+02 IFfM.GT. MSTOP) RETURN 

IF(CENTER(M) .LT.LTOL) GO TO 1+01 
I F (CENTER(M) . GT. HTOL) GO TO 301 
C WE REACH THIS POINT ONLY IF WE HAVE A HIT. 

C IF THE IDENTIFIED STANDARD MASS IS ALREADY IN THE FILE IGNORE 
C THIS FIND. 

DO 1+50 J=l , NUMSTD 

I F (CENTER (M) .EQ.STDDIS(J)) GO TO 800 
IF (STOMAS (J) .EQ. REFMAS (K) ) GO TO 800 
1+50 CONTINUE 

C IF THE NEW STANDARD MASS IS LARGER THAN THE LARGEST MASS ALREADY 
C IN THE FILE THEN WE CAN ADD IT TO THE END OF THE LIST. 

500 I F (CENTER (M) . LT. STDDI S (NUMSTD) ) GO TO 600 
NUMSTD=NUMSTD+1 
STDMAS (NUMSTD ) =RE FMAS ( K) 

STDD I S (NUMSTD )=C ENT ER(M) 

J=NUMSTD 
GO TO 800 

C THIS SECTION WILL INSERT AN IDENTIFIED STANDARD MASS WITHIN THE 
C STANDARD MASS LIST IN THE PROPER PLACE. THE MASSES MUST BE STORED IN 
C ASENDING ORDER. 

600 LIMIT=NUMSTD 

DO 625 J-l, LIMIT 

I F (CENTER (M) . GT. STDD IS (J) ) GO TO 625 

IT0P=1 60-J 

DO 650 LL=1 , ITOP 

JJ=160-LL 

JJJ=JJ+1 

STDDI S(JJJ)=STDDIS(JJ) 

650 STDMAS (JJJ)=STDMAS(JJ) 

STDD I S (J)=CENTER (M) 

STDMAS ( J ) =REFMAS ( K) 

NUMSTD=NUMSTD+1 
GO TO 800 
625 CONTINUE 

C THIS SECTION CHOOSES THE STANDARD MASSES THAT WILL BE USED TO 
C CALCULATE A NEW SET OF COEFFICIENTS. NOTE THAT VARIABLE J 
C POINTS TO WHERE IN THE STDD IS LIST THE LAST STANDARD MASS 
C WAS INSERTED OR ALREADY FOUND. 

800 IF (NUMSTD. EQ. 2) GO TO 1+01 
IF(STDDIS(J+l).NE.O.) J=J+1 
I F (J . LT.3) J-3 
DO 850 NEW=1 ,3 
INDEX=J-3+NEW 
X (NEW) “STDD l S ( INDEX) 

850 Y (NEW)=STDMAS (INDEX) 

NUM=3 
IDEG-2 
CALL POLYFT 
GO TO 301 
END 
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C JANUARY 29,1973 

C 

C THIS MODULE WILL ACCEPT CHANGES AND/OR DELETIONS TO THE STANDARD MASS LIST 
C 

C VARIABLES 

C 

C STDD I S ( 1 60 ) - THE DISTANCE VALUES OF PEAK CENTERS IDENTIFIED AS STANDARD MASSES 

C STDMAS (160) - THE SQUARE ROOT OF THE MASS VALUE OF IDENTIFIED STANDARD MASSES 

C NUMSTD - THE NUMBER OF IDENTIFIED STANDARD MASSES 

C REFMAS (160) - A TABLE OF THE MOST INTENSE LINES IN THE MARKER COMPOUND SPECTRUM 

C - THESE ARE THE MASSES SEARCHED FOR BY THE STANDARD MASS IDENTIFICATION 

C - ROUTINE 

C 

C xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE CHANGE 

DOUBLE PRECISION STDD I S ( 1 60) ,STDMAS(160) ,REFMAS(160) 

COMMON/PFKLST/ STDD I S , STOMAS , REFMAS , NUMSTD 

5000 WRITE(6,1) 

1 FORMAT ('0 '.'ENTER SUBSCRIPT') 

CALL INTGER(MM) 

IF(MM.EQ.O) go to 5001 
WRITE (6 ,3) 

3 FORMAT (IX, 'ENTER NEW MASS VALUE') 

C IF IT IS DESIRED TO DELETE A PEAK. ENTER ZERO FOR THE MASS VALUE 
READ (5,2) STDMAS(MM) 

2 FORMAT (F10.0) 

ST DMAS (MM) =DSQRT (STDMAS (MM) ) 

GO TO 5000 

C AFTER ALL CHANGES AND DELETIONS HAVE BEEN MADE, THE STANDARD MASS LIST MUST BE COMPRESSED. 

5001 M-0 

5002 M-M+l 

5003 I F (STDMAS (M) .NE.O.) GO TO 5002 
DO 5005 J«M, NUMSTD 

IF (STDMAS (J) .NE.O.) GO TO 5008 
5005 CONTINUE 
NUMSTD-M-1 
RETURN 

5008 DO 5010 J-M, NUMSTD 
STDD I S (J)-STDD I S (J+l ) 

5010 STDMAS (J)-STDMAS(J+1) 

GO TO 5003 
END 
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C JANUARY 29,1973 

C 

C THIS MODULE WILL DISPLAY AREAS OF THE CENTER LIST AND ADD ENTER I ES TO THE STANDARD 

C MASS LIST 

C 
C 
C 

C CENTER(1500) 

C 

C HEIGHT (1500) 

C NUMPK 

C STDD I S ( 1 60) 

C STDMAS ( 1 60) 

C NUMSTD 

C REFMAS(160) 

C 
C 

C TEMP 

C 

c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE ADD 

DOUBLE PRECISION STDD I S ( 1 60) ,STDMAS(l60) ,REFMAS(160) ,TEMP 
DOUBLE PRECISION CENTER 
COMMON/PFKLST/STDD I S , STDMAS , REFMAS , NUMSTD 
COMMON/PKLST /CENTER ( 1 500) ,HIGHT(1500) , NUMPK 
2 WRITE(6,30) 

30 FORMATCO' , 'ENTER DISTANCE AREA FOR DISPLAY') 

C ENTER ZERO TO RETURN TO THE CALLING MODULE. ENTER A VALUE LESS THAN 1. TO SKIP THE 
C DISPLAY SECTION AND GO DIRECTLY TO THE PEAK INPUT SECTION. 

READ (5,31 ) TEMP 

31 FORMAT (FI 0 .7) 

I F (TEMP . EQ. 0. ) GO TO 999 
I F (TEMP. LT. 1 . ) GO TO 41 
M=0 

1 M=M+ 1 

IF (CENTER (M) . EQ.O.) GO TO 1000 
IF(CENTER(M) .LT.TEMP) GO TO 1 
M=M-3 
K=M+6 

WRITE (6,32) 

C LIST THE RAW PEAK DATA IN THE VACINITY OF THE DESIRED PEAK 

32 FORMATCO' ,10X, 'DISTANCE' ,10X, 'RAW INTENSITY') 

IF(M.LT.l) M=1 

DO 40 J=M , K 

IF(CENTER(J).EQ.O.) GO TO 41 

40 WRITE(6,33) J, CENTER (j) ,HI GHT (J ) 

33 F0RMAT(4X, I3,4X,F8.4,12X,F7.1) 

41 WRITE(6,34) 

34 FORMATCO' , 'ENTER SUBSCRIPT') 

CALL INTGER(L) 

IF(L.EQ.O) GO TO 2 
WRITE (6,35) 

35 FORMAT (IX, 'ENTER MASS') 

READ (5, 36) TEMP 

36 FORMAT (FI 0. 0) 

TEMP=DSQRT (TEMP) 

NN=NUMSTD 

C PLACE THE ADDED PEAK IN THE PROPER PLACE IN THE STANDARD MASS LIST 
DO 7777 1-1 ,NN 

I F (TEMP . GT . STDMAS ( I ) ) TO TO 7777 

LTOP= 160-1 

DO 8000 J=1,LT0P 

MM=160-J 

MMM=MM+1 

STDD I S (MMM)=STDD I S (MM) 

8000 STDMAS (MMM)=STDMAS (MM) 

STDD I S (MM)=CENTER(L) 

STDMAS (MM) =TEMP 
NUMSTD=NUMSTD+1 
GO TO 2 
7777 CONTINUE 

C IF THE NEW PEAK IS LARGER THAN THE OTHER ENTERIES, ADD IT TO THE END OF THE LIST. 
NUMSTD=NUMSTD+1 
STDD I S(NUMSTD)=CENTER(L) 

STDMAS (NUMSTD) =TEMP 
GO TO 2 

1000 WRITE (6,37) CENTER (M-l ) 

37 FORMATCO' , 'THE LARGEST DISTANCE VALUE IN THE FILE I S' , F10.5) 

GO TO 2 

999 RETURN 
END 


VARIABLES 

- DISTANCE OF PEAK CENTERS FROM THE LEFT EDGE OF THE PHOTOPLATE 

- IN MILLIMETERS 

- MAXIMUM INTENSITY VALUE OF EACH PEAK 

- NUMBER OF PEAKS IN THE PEAK LIST TABLE 

- THE DISTANCE VALUES OF PEAK CENTERS IDENTIFIED AS STANDARD MASSES 

- THE SQUARE ROOT OF THE MASS VALUE OF IDENTIFIED STANDARD MASSES 

- THE NUMBER OF IDENTIFIED STANDARD MASSES 

- A TABLE OF THE MOST INTENSE LINES IN THE MARKER COMPOUND SPECTRUM. 

- THESE ARE THE MASSES SEARCHED FOR BY THE STANDARD MASS IDENTIFICATION 

- ROUTINE 

- TEMPORARY STORAGE 
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THIS MODULE WILL LIST THE STANDARD MASSES FOUND BY THE SEARCHING ROUTINE 


STDD I S ( 1 60) 
STDMAS (160) 
NUMSTD 
REFMAS(160) 


X (200) 

Y (200) 

N 

IDEG 

C(ll) 

LSTPAR ( 1 0) 
PFKMAS (160) 

STD (160) 

D I FF ( 1 60) 

PPM ( 1 6) 

CC(160, 1 1) 

INUM 

INDEX 
IHDEX 
I ORDER 

ISTRT 


M 


VARIABLES 

THE DISTANCE VALUES OF PEAK CENTERS IDENTIFIED AS STANDARD MASSES 
THE SQUARE ROOT OF THE MASS VALUE OF IDENTIFIED STANDARD MASSES 
THE NUMBER OF IDENTIFIED STANDARD MASSES 

A TABLE OF THE MOST INTENSE LINES IN THE MARKER COMPOUND SPECTRUM. 
THESE ARE THE MASSES SEARCHED FOR BY THE STANDARD MASS IDENTIFICATION 
ROUTINE 

UP TO 200 ABCISSA VALUES TO BE PASSED TO THE POLYNOMIAL FITTING 
ROUTINE 

UP TO 200 ORDINATE VALUES TO BE PASSED TO THE POLYNOMIAL FITTING 
ROUTINE 

NUMBER OF X AND Y VALUES PASSED 

DEGREE OF THE POLYNOMIAL (UP TO TENTH POSSIBLE) 

UP TO 11 COEFFICIENTS RETURNED BY THE POLYNOMIAL ROUTINE 
LIST OF PARAMETERS FROM THE PARLIST MODULE 

THE MASSES OF IDENTIFIED STANDARD MASSES CALCULATED BY USING THE 
SURROUNDING INUM STANDARD MASSES 

THE KNOWN EXACT MASS OF EACH IDENTIFIED STANDARD MASS 

THE DIFFERENCE IN MILLIMASS UNITS BETWEEN THE EXACT MASS AND THE 

CALCULATED MASS OF EACH IDENTIFIED STANDARD LINE 

THE DIFFERENCE IN PARTS PER MILLION BETWEEN THE EXACT MASS AND THE 

CALCULATED MASS OF EACH IDENTIFIED STANDARD LINE 

THE UP TO 160 SETS OF COEFFICIENTS CALCULATED BY THIS MODULE FOR 

EACH SET OF INUM STANDARD MASSES USED IN THE DISTANCE TO MASS 

CALCULATION 

THE NUMBER OF STANDARD MASSES USED FOR EACH SUCCESSIVE 

POLYNOMIAL FIT 

POINTER TO THE CC ARRAY 

THE TOTAL NUMBER OF SETS OF COEFFICIENTS 

THE NUMBER OF COEFFICIENTS RETURNED BY THE POLYNOMIAL FITTING 
ROUTINE 

POINTER TO THE LOWEST OF THE INUM VALUES TO BE TRANSFERED TO 
THE POLYNOMIAL FITTING ROUTINE 

POINTER TO THE HIGHEST OF THE INUM VALUES TO BE TRANSFERED TO 

THE POLYNOMIAL FITTING ROUTINE 

THE STDD IS AND ST DMAS ARRAY POINTER 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE STDLST 

DOUBLE PRECISION X(200) ,Y(200) ,C(1 1 ) 

DOUBLE PRECISION STDMAS ( 1 60) , STDD I S ( 1 60) , PPM ( 1 60) , REFMAS ( 1 60) 

DOUBLE PRECISION PFKMAS ( 1 60) , STD ( 1 60) ,D I FF( 160) ,CC ( 160, 1 1 ) 

COMMON/PFKLST /STDD I S , STDMAS , REFMAS , NUMSTD 
COMMON/F IT/X, Y,N, IDEG.C 
COMMON/OUTLST/ PFKMAS , STD , D I FF , PPM , CC 
COMMON/PAR/LSTPAR (10) 

WRITE (6, 11) 

SET THE NUMBER OF STANDARD MASSES USED, AND THE DEGREE OF THE CONVERSION POLYNOMIAL 
INUM=LSTPAR(7) 

LIMDEG«LSTPAR(8) 

CLEAR VARIABLES 

DO 100 M=1 ,160 
PFKMAS (M)*0. 

STD(M)-0. 

PPM (M)*0. 

D I FF (M)“0. 

DO 100 J-1,11 
C(J)«0. 

100 CC(M,J)-0. 

INUM IS THE NUMBER OF POINTS USED BY THE DISTANCE TO MASS CONVERSION POLYNOMIAL 
JNUM IS ONE LESS THAN INUM, AND KNUM IS ONE HALF INUM. 

JNUM-INUM-1 

KNUM-INUM/2 

THIS SECTION IS FOR LESS THAN INUM IDENTIFIED STANDARD MASSES 
I F (NUMSTD. GT. INUM) GO TO 1000 

IF THERE ARE LESS THAN INUM IDENTIFIED STANDARD MASSES, SEND THEM ALL TO CALCULATE 
A POLYNOMIAL OF DEGREE NUMSTD-2'; PROVIDED THAT THE DEGREE DOES NOT EXCEED LIMDEG 
N=NUMSTD 

DO 500 M-l, NUMSTD 
X(M)«STDDIS(M) 

500 Y (m)=stdmas (m) 

IDEG-NUMSTD-2 

IF(IDEG.GT. LIMDEG) IDEG-LIMDEG 
CALL POLYFT 
I ORDER** I DEG+ 1 

STORE THE RETURNED SET OF COEFFICIENTS IN THE CC ARRAY, IN THIS CASE THERE WILL BE 
ONLY ONE SET. 

DO 501 M*1 , IORDER 
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C JANUARY 26,1973 
C 

501 CC (1 ,M)=C (M) 

I HDEX-1 
M=1 

GO TO 3000 

C THIS SECTION IS FOR MORE THAN INUM IDENTIFIED STANDARD MASSES 

1000 INDEX-1 

I START* 1 

1001 ITOP- I START+JNUM 

I F(LTOP.GT.NUMSTD) GO TO 2000 
MM=0 

C ALL THE STANDARD MASSES ARE SENT TO THE POLYNOMIAL FITTING ROUTINE IN GROUPS OF INUM. 

C AFTER ALL THE MASSES HAVE BEEN USED FOR POLYNOMIAL CALCULATION IN THIS MANNER, THERE 
C WILL BE NUMSTD-JNUM SETS OF COEFFICIENTS STORED IN THE CC ARRAY. THESE SETS OF COEFFICIENTS 
C REPRESENT THE BEST POLYNOMIAL FITS FOR DISTANCE TO MASS CONVERSIONS IN THE NUMSTD-JNUM 
C OVERLAPPING AREAS OF THE STANDARD MASS LIST. 

DO 300 M-l START, I TOP 
MM=MM+ 1 

X (MM)-STDD I S (M) 

300 Y(MM)=STDMAS(M) 

N- I NUM 
I DEG-L IMDEG 

C STORE COEFFICIENTS IN THE PROPER PLACE IN THE CC ARRAY 
CALL POLYFT 
DO 350 M-l ,11 
350 CC ( I NDEX,M)-C (M) 

I NDEX- I NDEX+1 
I START®! START+1 
GO TO 1001 

C AFTER ALL SETS OF COEFFICIENTS HAVE BEEN CALCULATED, USE THEM TO DETERMINE THE ERROR 
C IN THE DISTANCE TO MASS CONVERSION. LIST THE DIFFERENCE BETWEEN THE KNOWN MASSES OF 
C THE REFERENCE PEAKS AND THE MASSES CALCULATED FOR THESE REFERENCE LINES USING THE SAME 
C METHOD OF CALCULATION AS FOR UNKNOWN LINES. 

2000 I NDEX-1 

IHDEX-NUMSTD-JNUM 

M-l 

3000 IF(M.GT.NUMSTD) GO TO 5000 
IF(M.LT.KNUM) I NDEX*=1 
I F ( INDEX.GT. IHDEX) INDEX-IHDEX 
PFKMAS (M) *»CC ( INDEX, 1 ) 

DO 4000 K- 1 ,10 
KK-K+1 

4000 PFKMAS (M)-PFKMAS (M)+CC (INDEX, KK) *STDD I S (M) **K 
M-M+l 

INDEX-1 NDEX+1 
GO TO 3000 

5000 DO 4500 K-l.NUMSTD 

PFKMAS (K) -PFKMAS (K)*PFKMAS (K) 

STD(K)=STDMAS(K)*STDMAS(K) 

D I FF (K)= (PFKMAS (K) -STD (K) )*1 000. 

4 500 PPM ( K) -DABS ( (D I FF ( K) /STD ( K) ) * 1 000 . ) 

C THIS IS THE LISTING SECTION 
WRITE (6,201 ) 

WRITE (6 , 202) 

WRITE (6,203) 

DO 4001 M-l ,NUMSTD 

4001 WRITE (6 ,204) M,STDDIS(M) ,PFKMAS(M) ,STD(M) ,DIFFCM) ,PPM(M) 

RETURN 

1 FORMAT ('1') 

201 FORMAT (' O' ,1 OX, 'D I STANCE' , 1 OX, 'CALCULATED' , 1 OX, 'EXACT MASS' , 10X, 'DIFFERENCE' , 10X, 'DIFFERENCE' ) 

202 FORMAT ( ' ' ,31X, 'MASS' ,35X, ' IN MMU ' , 1 4X , ' I N PPM') 

203 FORMAT (' ') 

204 F0RMAT(3X, I3,5X,F8.4,9X,F10.5,10X,F10.5J3X,F5.2,15X,F5.2) 

END 
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C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


JANUARY 30,1973 

THIS MODULE ACCEPTS ADDITIONS TO THE ELEMENT LIST USED FOR FITTING COMPOSITIONS 
TO THE UNKNOWN PEAKS. IT ALSO WILL ACCEPT A CHANGE IN THE ELEMENT FITTING TOLERANCE 
SUPPLIED BY THE PARAMETER SETTING MODULE. IT WILL, AT THE USERS OPTION, WRITE ALL THE DATA 
OBTAINED SO FAR INTO A DISK FILE OR WILL CALL THE MODULES TO DO THE ELEMENT FITTING FOR 
LISTING ON THE CONVERSATIONAL TERMINAL. 


VARIABLES 


C ENTER ( 1 500) 

H I GHT ( 1 500) 
NUMPK 

STDD I S ( 1 60) 
STDMAS(160) 
NUMSTD 
REFMAS(160) 


LSTPAR ( 1 0) 
PFKMAS (160) 

STD (160) 

D I FF ( 1 60) 

PPM (160) 

CC (160, 1 1 ) 

MASS (12) 
NAME ( 12) 
NUMBER (12) 
NTOL 
NUM 


- DISTANCE OF PEAK CENTERS FROM THE LEFT EDGE OF THE PHOTOPLATE 

- IN MILLIMETERS 

- MAXIMUM INTENSITY VALUE OF EACH PEAK 

- NUMBER OF PEAKS IN THE PEAK LIST TABLE 

- THE DISTANCE VALUES OF PEAK CENTERS IDENTIFIED AS STANDARD MASSES 

- THE SQUARE ROOT OF THE MASS VALUE OF IDENTIFIED STANDARD MASSES 

- THE NUMBER OF IDENTIFIED STANDARD MASSES 

- A TABLE OF THE MOST INTENSE LINES IN THE MARKER COMPOUND SPECTRUM. 

- THESE ARE THE MASSES SEARCHED FOR BY THE STANDARD MASS IDENTIFICATION 

- ROUTINE 

- LIST OF PARAMETERS FROM THE PARLIST MODULE 

- THE MASSES OF IDENTIFIED STANDARD MASSES CALCULATED BY USING THE 

- SURROUNDING STANDARD MASSES 

- THE KNOWN EXACT MASS OF EACH IDENTIFIED STANDARD MASS 

- THE DIFFERENCE IN MILLIMASS UNITS BETWEEN THE EXACT MASS AND THE 

- CALCULATED MASS OF EACH IDENTIFIED STANDARD LINE 

- THE DIFFERENCE IN PARTS PER MILLION BETWEEN THE EXACT MASS AND THE 

- CALCULATED MASS OF EACH IDENTIFIED STANDARD LINE 

- THE UP TO 160 SETS OF COEFFICIENTS CALCULATED FOR EACH SET OF 

- STANDARD MASSES USED IN THE DISTANCE TO MASS CALCULATION 

- UP TO TWELVE MASSES OF ELEMENTS USED FOR FITTING 

- A FOUR CHARACTER NAME FOR UP TO TWELVE ELEMENTS 

- THE MAXIMUM NUMBER OF EACH ELEMENT USED IN ELEMENT FITTING 

- THE TOLERANCE IN MILLIMASS UNITS USED FOR ELEMENT FITTING 

- THE NUMBER OF ELEMENTS IN THE MASS ARRAY 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE ELELST 

DOUBLE PRECISION CENTER, STDMAS(160) ,STDD I S (1 60) , PPM ( 1 60) ,REFMAS(l60) 

DOUBLE PRECISION PFKMAS(l60) ,STD(160) ,DIFF(160) ,CC(160, 1 1 ) , MASS ,MAS ( 1 2 ) 

INTEGER NAM (12) , NUMB (12) 

COMMON/PFKLST/STDD I S , STDMAS , REFMAS , NUMSTD 
COMMON/PAR/ LSTPAR (10) 

COMMON/OUTLST/ PFKMAS , STD ,D I FF, PPM ,CC 
COMMON/PKLST/ CENTER(1500) ,HE I GHT(1 500) , NUMPK 
COMMON/ELEM/ MASS(12) ,NAME(12) ,NUMBER(12) , NTOL, NUM 
DATA MAS/1. 007821*6,12. ,13- 0033554, 1^.0030732, 15. 9949141 ,7*0/ 

DATA NAM/' H 1 , ' C 1 2 ' , ' C13 ','N14 ','016 ',7*0/ 

DATA NUMB/2*O,2,3,4,7*0/ 

DATA I YES, IBLK/'Y' , 1 '/ 

DO 50 M* 1,12 
MASS (M)=MAS (M) 

NAME (M)=NAM (M) 

50 NUMBER (M)=NUMB (M) 

1-5 


NUM=5 
NTOL-5 
WRITE (6, 1 00) 

READ (5,51 ) J 
I F ( J . NE . I YES) GO TO 2001 

C ENTER NEW VALUE FOR ELEMENT FITTING TOLERANCE OR ENTER ZERO TO RETAIN OLD VALUE 
WR ITE (6, 1 06) NTOL 
CALL INTGER(II) 

I F ( I I .NE.O) NTOL® I I 

C ENTER ADDITIONAL ELEMENTS USED FOR FITTING 
1 1*1+1 

I F( I .EQ. 13) GO TO 2000 
WR ITE (6, 1 02) 

READ (5,52) NAME ( I ) 

I F (NAME ( I ) . EQ. IBLK) GO TO 1999 
WRITE (6, 104) 

READ (5, 53) MASS ( I ) 

WRITE (6, 105) 

CALL INTGER(NUMBERO)) 

GO TO 1 

1999 NAME ( I )=0. 

2000 NUM= I - 1 

2001 WRITE(6,103) 

READ (5,50 I FLAG 

I F ( I FLAG . EQ.. I YES) GO TO 3000 

C WRITE ALL DATA INTO DISK FILE FOR ACCESS BY THE BATCH ELEMENT FITTING ROUTINE 

WR I TE (3) NUMPK, (CENTER ( I ) , H I GHT ( I ) , I = 1 , NUMPK) , STDD I S , STDMAS , REFMAS , NUMSTD , • 
XLSTPAR , PFKMAS , STD , D I FF, PPM ,CC , MASS , NAME , NUMBER , NTOL , NUM 
STOP 
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C JANUARY 30,1973 
C 

C CALL THE ELEMENT FITTING ROUTINE FOR CALCULATION AND DISPLAY ON THE CONVERSATIONAL 
C TERMINAL 

3000 CALL ELEMFT 
RETURN 

100 FORMAT ( 1 0 1 , 1 DO YOU WISH TO REVIEW ELEMENT FITTING PARAMETERS? 1 ) 

102 FORMAT (' O' , 'ENTER ISOTOPE NAME') 

103 FORMAT ( ' O' , ' I S THIS TERMINAL OUTPUT ?') 

10A FORMAT (' ',' ENTER MASS') 

105 FORMAT (' ',' NUMBER FOR THIS ISOTOPE') 

106 F0RMAT(1X, 'THE TOLERANCE FOR ELEMENT FITTING IS', 13,' MMU' ) 

51 FORMAT (Al) 

52 FORMAT (AA) 

53 FORMAT (F10. 0) 


END 
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C 

C 

C 

C 
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THIS MODULE WILL CALCULATE AN EXACT MASS FOR EACH UNKNOWN PEAK IN THE PEAK LIST. IT 
WILL ALSO DETERMINE ALL ELEMENTAL COMPOSITIONS THAT ARE WITHIN A SPECIFIED TOLERANCE 
OF THE EXACT MASS FROM THE GIVEN LIST OF ELEMENTS. 


VARIABLES 


CENTER(1500) 

HIGHT (1500) 
NUMPK 

STDD I S ( 1 60) 
STDMAS(160) 
NUMSTD 
REFMAS ( 1 60) 


LSTPAR(IO) 
PEKMAS (160) 

STD (160) 

D I FF ( 1 60) 

PPM( 16) 

CC ( 1 60 ,11) 

MASS (12) 

NAME (12) 
NUMBER( 1 2) 
NTOL 
NUM 

CALMAS (391 ) 

CLCMAS ( 1 500) 
ATNUM ( 1 2) 

COUNT 1 
C0UNT2 


- DISTANCE OF PEAK CENTERS FROM THE LEFT EDGE OF THE PHOTOPLATE 

- IN MILLIMETERS 

- MAXIMUM INTENSITY VALUE OF EACH PEAK 

- NUMBER OF PEAKS IN THE PEAK LIST TABLE 

- THE DISTANCE VALUES OF PEAK CENTERS IDENTIFIED AS STANDARD MASSES 

- THE SQUARE ROOT OF THE MASS VALUE OF IDENTIFIED STANDARD MASSES 

- THE NUMBER OF IDENTIFIED STANDARD MASSES 

- A TABLE OF THE MOST INTENSE LINES IN THE MARKER COMPOUND SPECTRUM. 

- THESE ARE THE MASSES SEARCHED FOR BY THE STANDARD MASS IDENTIFICATION 

- ROUTINE 

- LIST OF PARAMETERS FROM THE PARLIST MODULE 

- THE MASSES OF IDENTIFIED STANDARD MASSES CALCULATED BY USING THE 

- SURROUNDING STANDARD MASSES 

- THE KNOWN EXACT MASS OF EACH IDENTIFIED STANDARD MASS 

- THE DIFFERENCE IN MILLIMASS UNITS BETWEEN THE EXACT MASS AND THE 

- CALCULATED MASS OF EACH IDENTIFIED STANDARD LINE 

- THE DIFFERENCE IN PARTS PER MILLION BETWEEN THE EXACT MASS AND THE 

- CALCULATED MASS OF EACH IDENTIFIED STANDARD LINE 

- THE UP TO 160 SETS OF COEFFICIENTS CALCULATED FOR EACH SET OF 

- STANDARD MASSES USED IN THE DISTANCE TO MASS CALCULATION 

- UP TO TWELVE MASSES OF ELEMENTS 

- A FOUR CHARACTER NAME FOR UP TO TWELVE ELEMENTS 

- THE MAXIMUM NUMBER OF EACH ELEMENT USED IN ELEMENT FITTING 

- THE TOLERANCE IN MILLIMASS UNITS USED FOR ELEMENT FITTING 

- THE NUMBER OF ELEMENTS IN THE MASS ARRAY 

- A LIST OF THE MASSES OF THE PEAKS IN THE MARKER COMPOUND THAT 

- WILL BE REMOVED FROM THE CLCMAS TABLE IF FOUND 

- A CALCULATED MASS FOR EACH PEAK IN ARRAY CENTER 

- THE NUMBER OF ATOMS OF EACH ELEMENT THAT FITS A PARTICULAR 

- ELEMENTAL COMPOSITION 

- A COUNTER THAT CONTROLS THE OUTPUT LISTING 

- THE NUMBER OF ELEMENTAL COMPOSITIONS FOUND SO FAR FOR THE 

- THE PARTICULAR CLCMAS VALUE UNDER CONSIDERATION 


XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 


SUBROUTINE ELEMFT 

DOUBLE PRECISION CENTER, STDMAS (1 60) , STDD I S (1 60) ,PPM(160) ,REFMAS(l60) 

DOUBLE PRECISION PFKMAS ( 1 60) , STD ( 160) ,D I FF(1 60) ,CC ( 1 60, 1 1 ) 

DOUBLE PRECISION CLCMAS (1 500) ,D IF, HTOL.LTOL.CALMAS (391 ) .MASS 
INTEGER C0UNT1 .C0UNT2, ATNUM (12) 

REAL MAX 

COMMON/PFKLST /STDD I S , STDMAS , REFMAS , NUMSTD 
COMMON/PAR/ LSTPAR(IO) 

COMMON/OUTLST/ PFKMAS , STD , D I FF , PPM ,CC 
COMMON/PKLST/ CENTER(1500) ,HIGHT(1500) , NUMPK 
COMMON/ELEM/ MASS (12) .NAME (12) .NUMBER (12) , NTOL, NUM 
COMMON/COMP/CLCMAS , ATNUM, COUNT 1 ,C0UNT2,DIF, IMAS 
C CLEAR CLCMAS ARRAY 
DO 50 1-1,1500 
50 CLCMAS ( I )=0. 

C SET THE NUMBER OF STANDARD MASSES USED FOR POLNOMIAL FITTING IN MODULE STDLST 
I NUM-LSTPAR (7) 

C CALL THE MARKER COMPOUND LISTING MODULE 
CALL STDLST 

C INITALIZE THE COUNTER THAT CONTROLS THE OUTPUT FORMAT 
C0UNT1-A5 

C CALCULATE THE EXACT MASS OF EACH PEAK USING THE PROPER SET OF POLYNOMIAL COEFFICIENTS 
C FROM MODULE STDLST. INDEX IS THE CC ARAY POINTER, IHDEX IS THE POINTER TO THE LAST 
C SET OF COEFFICIENTS IN THE CC ARRAY. 

1012 INDEX-1 

I HDEX-NUMSTD- I NUM+1 
I F( IHDEX.LT. 1) IHDEX-1 
M=1 

1 DO 1000 J— 1 , NUMSTD 

I F (STDD I S (J) . LT. CENTER (M) ) GO TO 1000 
GO TO 1001 

1000 CONTINUE 

1001 INDEX«J-( INUM/2) 

I F (J . LE . I NUM/2) INDEX-1 
IF(INDEX.GT. IHDEX) INDEX- IHDEX 
CLCMAS (M)-CC( I NDEX.l) 

DO 1005 K— 1,10 
KK-K+1 

1 005 CLCMAS (M) -CLCMAS (M) +CC (INDEX, KK) *CENTER (M) **K 


C FEBRUARY !, 1973 MICHAEL ROMIEZ 

C 

M=M+I 

IF(M.GT.NUMPK) GO TO 2000 
GO TO 1 

2000 DO 1010 M=1 , NUMPK 

1010 CLCMAS (M)=CLCMAS (M) *CLCMAS (M) 

C READ IN MASSES TO BE STRIPED OUT 
READ (A) CALMAS 

C STRIP OUT BACKGROUND AND MARKER COMPOUND PEAKS 
TOL=FLOAT(NTOL)*.0Ol 
J=0 
10 J=J+1 

IF(J.GT.NUMKP) GO TO 2500 
DO 2001 M=1 ,391 
HTOL=CALMAS(M)+TOL 
LTOL=CALMAS (M) — TOL 
I F (CLCMAS (J) . LT. LTOL) GO TO 10 
I F (CLCMAS (J) . GT. HTOL) GO TO 2001 
CLCMAS (J)=0. 

GO TO 10 

2001 CONTINUE 
GO TO 10 

C COMPRESS CLCMAS ARRAY 
2500 M=0 

100 M=M+1 

170 I F (CLCMAS (M) . NE. 0. ) GO TO 100 
DO 150 J=M, NUMPK 
CLCMAS (J)=CLCMAS (J+l ) 

CENTER(J)=CENTER (J+l ) 

150 H I GHT ( J ) =H I GHT (J+l ) 

DO 157 K=M, NUMPK 
I F (CLCMAS (K) . EQ.O. ) GO TO 167 
GO TO 170 

167 CONTINUE 

C IF THERE ARE NO PEAKS IN THE MASS SPECTRUM, RETURN TO THE CALLING MODULE 
IF(CLCMAS(1) .EQ.O) RETURN 
C CALCULATE NEW NUMPK VALUE 
2700 DO 2800 L=1 .NUMPK 

I F (CLCMAS (L).NE.O.) GO TO 2800 
LL=L- 1 
TO TO 2900 
2800 CONTINUE 
2900 NUMPK=LL 
C NORMALIZE INTENSITY VALUES 
MAX=H I GHT ( 1 ) 

00 6950 J=2, NUMPK 
I F (H I GHT (J ) . GT . MAX ) MAX=HIGHT(J) 

4950 CONTINUE 

I F (MAX. EO. 0) GO TO 3005 
DO 4951 J-l .NUMPK 

4951 HIGHT(J)=(HIGHT(J) /MAX ) * 1 00 . 

C CALCULATE ELEMENTAL COMPOSITIONS 

3005 NUM3=NUMBER(3) 

NUM4=NUMBER(4) 

NUM5=NUMBER(5) 

NUM6=NUMBER(6)+1 

NUM7“NUMBER(7)+1 

num8=number(8)+i 

NUM9=NUMBER (9)+l 

NUM10=NUMBER(10)+1 

NUM1 1=NUMBER(1 1 )+l 

NUM12=NUMBER(12)+1 

DO 4000 M=l, NUMPK 

HTOL=CLCMAS(M)+TOL 

LTOL=CLCMAS(M)-TOL 

C0UNT2=0 

NUM2=I D INT (HTOL/1 2. )+l 
NUMl=NUM2*2+8 
00 4050 112=1 ,NUM12 
1112=112-1 
ATNUM ( 1 2 ) = I I 12 
DO 4050 1 1 1=1 ,NUM1 1 
1111=111-1 
ATNUM(1 l)=llll 
DO 4050 1 10=1 ,NUM10 
1110=110-1 
ATNUM (1 0 ) = I 1 10 
DO 4050 19-1 ,NUM9 
1 19=19-1 
ATNUM(9)=I 19 
DO 4050 18=1, NUM8 
1 18=18-1 
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C FEBRUARY I, 1973 MICHAEL ROMIEZ 

C 

ATNUM(8)=I 18 
DO 4050 I7-1.NUM7 
1 17-17-1 
ATNUM(7)=I 17 
DO 4050 16=1, NUM6 
1 16=16-1 
ATNUH (6 ) = I 16 
DO 4050 I5-1.NUM5 
115=15-1 
ATNUM(5)=I 15 
DO 4050 14=1, NUM4 
1 14=14-1 
ATNUM(4)=I 14 
DO 4050 13=1 ,NUM3 
113=13-1 
ATNUM(3) = I 13 
DO 4050 12=1 ,NUM2 
I 12=12-1 
AT N UM ( 2 ) = I 12 
SUM-0. 

DO 4100 I 1=2, NUM 
4100 SUM-SUM+ATNUM (II) *MASS (II) 

C ITEST IS AN INTEGER VALUE WHICH IS USED TO FIT THE ELEMENT HYDROGEN INTO THE 
C ELEMENTAL COMPOSITIONS. IF ITEST IS LESS THAN ZERO THIS MEANS THAT SUM 
C IS GREATER THAN CLCMAS(M) PLUS THE TOLERANCE AND WE WILL DISMISS THIS COMPOSITION. 

C WHEN ITEST IS ZERO OR POSITIVE WE HAVE A POSSIBLE COMPOSITION AND WE MAKE FURTHER 
C TESTS 

I T E ST= I D I NT (HTOL-SUM) 

I F ( ITEST.GT. NUM1 ) GO TO 4050 
I F ( ITEST) 4050,4051,4051 

4051 ATNUM(1 ) = I T E ST 
SUM-SUM+MASS ( 1 ) *ATNUM ( 1 ) 

GO TO 4053 

C WE CONTINUE TO ADD HYDROGEN ATOMS TO THE COMBINATIONS OF OTHER ELEMENTS UNTIL WE HAVE 
C A FIT OR WE HAVE TOO HIGH OF A TOTAL MASS FOR A FIT. 

4052 SUM-SUM+MASS ( I ) 

ATNUM(1 )=ATNUM (1 )+l 

4053 I F (SUM. GT. HTOL) GO TO 4050 
IF SUM.LT.LTOL) GO TO 4052 
C0UNT2-C0UNT2+1 

0 1 F=(CLCMAS (M) -SUM) *1000. 

IMAS-M 
CALL OUTPUT 
4050 CONTINUE 

C IF THERE WERE NO ELEMENTAL COMPOSITIONS FOUND FOR THIS MASS, LIST THIS INFORMATION 
C AS PART OF THE OUTPUT. 

IMAS-M 

I F (C0UNT2. EQ. 0) CALL OUTPUT 
4000 CONTINUE 
RETURN 
END 
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C FEBRUARY 1,1973 MICHAEL ROMIEZ 

C 

C THIS MODULE LISTS ELEMENTAL COMPOSITIONS FOUND BY THE ELEMENT FITTING ROUTINE 
C 

C VARIABLES 

C 

C CENTER (1500) - DISTANCE OF PEAK CENTERS FROM THE LEFT EDGE OF THE PHOTOPLATE 

C - IN MILLIMETERS 

C HIGHT (1500) - MAXIMUM INTENSITY VALUE OF EACH PEAK 

C NUMPK - NUMBER OF PEAKS IN THE PEAK LIST TABLE 

C MASS (12) - UP TO TWELVE MASSES OF ELEMENTS 

C NAME (12) - A FOUR CHARACTER NAME FOR UP TO TWELVE ELEMENTS 

C NUMBER(12) - THE MAXIMUM NUMBER OF EACH ELEMENT USED IN ELEMENT FITTING 

C NTOL - THE TOLERANCE IN MILLIMASS UNITS USED FOR ELEMENT FITTING 

C NUM - THE NUMBER OF ELEMENTS IN THE MASS ARRAY 

C ATNUM(12) - THE NUMBER OF ATOMS OF EACH ELEMENT THAT FITS A PARTICULAR 

C - ELEMENTAL COMPOSITION 

C COUNT 1 - A COUNTER THAT CONTROLS THE OUTPUT LISTING 

C C0UNT2 - THE NUMBER OF ELEMENTAL COMPOSITIONS FOUND SO FAR FOR THE 

C - THE PARTICULAR CLCMAS VALUE UNDER CONSIDERATION 

C 

c xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

SUBROUTINE OUTPUT 

DOUBLE PRECISION CENTER, MASS, CLCMAS(1500) ,D IFF 
INTEGER ATNUM (12), COUNT 1 ,C0UNT2 
COMMON/ELEM/MASS ( 12) ,NAME(12) ,NUMBER(12) , NTOL, NUM 
COMMON/PKLST/ CENTER(ISOO) ,HIGHT(1500) .NUMPK 
COMMON/COMP/CLCMAS, ATNUM, C0UNT1 .C0UNT2 ,D I FF, IMAS 
C TEST NEED FOR A NEW PAGE HEADING 
I F (COUNT 1 .GE.^5) GO TO 1002 

C IF C0UNT2 IS ZERO, THERE HAVE BEEN NO ELEMENTAL FITS FOR THIS MASS 
500 I F (C0UNT2 . EQ.. 0) GO TO 1001 
C IF THIS IS A NEW MASS WE SKIP A SPACE 
IF(C0UNT2.EQ. 1) GO TO 1003 

1000 WRITE (6,303) CENTER ( IMAS) ,HIGHT ( IMAS) , CLCMAS (IMAS) ,DIFF, (ATNUM(M) ,M“1 ,NUM) 

COUNT 1 “COUNT 1 + 1 

RETURN 

1001 WR ITE (6,30*0 C ENTER (I MAS) ,HIGHT( I MAS) .CLCMAS (IMAS) ,NTOL 
COUNT 1 “COUNT 1+2 

RETURN 

1002 WRITE (6, 100) 

WRITE(6,300) 

WRITE (6,301 ) (NAME (M) ,M“1 ,NUM) 

WRITE (6,302) 

COUNT 1-5 
GO TO 500 

1003 WRITE(6,203) 

COUNT 1 “COUNT 1 + 1 
GO TO 1000 

203 FORMAT ( ' ' ) 

1 00 FORMAT ( ' 1 ' ) 

300 FORMAT (6X, 'DISTANCE' ,6X, 'RELATIVE' ,7X, 'MASS' ,7X, 'DIFFERENCE') 

301 F0RMAT(20X, ' INTENSITY' ,19X,' IN MMU' ,5X, 12A5) 

302 FORMAT ( ' 0 ' ) 

303 FORMAT (6X,F8. 4, 7X,F5. 1 , 5X , F 1 0 . 5 , 8X , F5 - 2 , 3X , 1 2 1 5 ) 

304 FORMAT ( ' 0 ' , 5X , F8 . 4 , 7X , F5 . 1 , 5X , F 1 0 . 5 , 1 0X , ' NO ELEMENTAL COMPOSITIONS FOUND WITHIN', 13,' MMU') 
END 
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